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ABSTRACT 

The structures of globular clusters (GCs) reflect their dynamical states and 
past histories. High-resolution imaging allows the exploration of morphologies 
of clusters in other galaxies. Surface brightness profiles from new Hubble Space 
Telescope observations of 34 globular clusters in M31 are presented, together 
with fits of several different structural models to each cluster. M31 clusters 
appear to be adequately fit by standard King models, and do not obviously 
require alternate descriptions with relatively stronger halos, such as are needed 
to fit many GCs in other nearby galaxies. The derived structural parameters 
are combined with corrected versions of those measured in an earlier survey to 
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construct a comprehensive catalog of structural and dynamical parameters for 
M31 GCs with a sample size similar to that for the Milky Way. Clusters in M31, 
the Milky Way, Magellanic Clouds, Fornax dwarf spheroidal and NGC 5128 define 
a very tight fundamental plane with identical slopes. The combined evidence 
for these widely different galaxies strongly reinforces the view that old globular 
clusters have near-universal structural properties regardless of host environment. 

Subject headings: galaxies: star clusters — galaxies: individual (M31) — globular 
clusters: general 



1. Introduction 

Globular clusters (GCs) have long been recognized as unique dynamical laboratories and 
unique markers of galaxy history. The integrated properties of GCs, such a s age and metallic- 
ity, a re believed to reflect conditions in the early stages of galaxy formation (IBrodie Strader 
20061 ). The spatial structure and kinematics of GCs reflect both their formation conditions 
and dynamical evolution within the tidal fields of their host galaxies. 

Structures of GCs as described by surface brightness or velocity profiles can be fit by 
a number of different models. These models can be used to compute numerous parameters 
describing various aspects of cluster structure, but four independent parameters suffice to 
completely define a single-mass model. It is helpful to think of cluster structures as being 
defined by spatial scale, spatial concentration, total luminosity, and a relation be t ween lumi- 



nosity and mass, b ut many different parameter choices are possible. iDjorgovskil (119951 ) and 



McLaughlin! (120001 ) used formally different, but entirely equivalent, sets of structural parame- 



ters to show that globular clusters in the Milky Way do not occupy the full four- dimensional 



appears to be true for clusters in M 31 (IDiorgovs 



Barmbv et al.l l2002h . in NGC 5128 ([Harris et al 



20021 ) 



Dubat 


i & Grillmair 


1997 


M33 ( 


Larsen et al. 


2002) 



This evidence indicates that the formation and evol ution process es for normal globular clus- 
ters were similar in a wide variety of environments (jHarrisI 120031 ) . 



1 Based on observations made with the NASA/ESA Hubble Space Telescope, obtained at the Space 
Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, 
Inc., under NASA contract NAS 5-26555. These observations are associated with programs #8664, #9719, 
# 9767 and #10260. 
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The globular cluster fundamental plane is simil ar in character, but different in detail, 
to that of elliptical galaxies (for an introduction, see lBurstein et al.lll997l ). There have been 
several suggestions that at least some of the most m assive GCs in the Milky Way and M31 



Tsuchiya et al 



2003: Hilker fc Richtler 



are in fact the nuclei of destroyed dwarf galaxies (e.g._ 
20001 ). or that they are related to ultra-compact dwarf galaxies ( jHasegan et alil2005h . Both 
dwarf galaxy nuclei and UCDs occupy somewhat different regions of the fundamental plane 
than do GCs, but structural properties have been measured for only a few objects in this 
transition mass region between globulars and galaxies. Combining samples of clusters from 
a number of nearby galaxies — where structures can be accurately measured with high- 
resolution imaging from the Hubble Space Telescope (HST) — is a way to improve the sample 
statistics particularly for the upper end of the GC mass range. 

To further understand the globular cluster fundamental plane, we have carried out a 
survey of m assive globulars i n M31 and NGC 5128 using the HST Advanced Camera for Sur- 



veys fA CS; iFord et a 



2002). Other papers from this program discuss the NGC 5128 cluster 
the tec hniques for measurement of those clusters' structural pa- 



sample (IHarris et al.L , ...... 

rameters ^McLaughlin et al. 2007a ). and the compari son between massive glob ulars, nuclear 



star clusters, and 'dwarf- globular transition objects' (jMcLaughlin et al.ll2007bl ). This paper 
presents the new M31 observations and structural parameters, compiles a publicly- available 
catalog of measurements of structural parameters for the largest-ever sample of M31 clusters, 
and compares the resulting 'fundamental planes' of globular cluster parameters in different 
galaxies. Throughout this paper we assume a distance to M31 of 780 kpc (1 pc subtends 



0'/26) and use the cluster nomenclature of the Revised Bologna Catalogue (IGalleti et al. 
20061 1 . 



2. Observational material 

We have used material from several different HST programs. Program GO- 10260 (PI 
W. Harris) obtained ACS Wide Field Channel (WFC) images centered on six clusters in M31 
and twelve clusters in NGC 5128, selected to be those having My < —10.5 and (for M31) 
no previous HST observations. Exposure times were 2370 s (one orbit split into three dither 
positions) in both the F606W and F814W filters. Most of the GO-10260 M31 observations 
were made in the period 2004 September 29-2004 October 01. Program GO-9719 (PI T. 
Bridges) obtained ACS High Resolution Channel (HRC) images of three M31 clusters using 
exposure times of 2020 s in F606W and 2680 s in F814W in 2003 August and September. 
Program GO-8664 (PI W. Harris) tar geted 24 M31 clus ters with a range of properties (as 



well as NGC 5128 clusters, discussed in IHarris et al.ll2002l ) for imaging with STIS in snapshot 
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mode; exposure times were 480 s in the unfiltered 50CCD/CL imaging mode. These data 
were obtained in fall 2000 and winter 2001. Of the 24 STIS targets, one (B329) is clearly a 
galaxy, and two more (B391 and B340) had faulty coordinates and were not on the images. 
Cluster B396 was bisected by the edge of the STIS image but the data were still usable. 

The cluster G001 (Mayall II) is one of the most massive in M31 and has been the 
subject of numerous recent investigations. It is often likened to the Milky Way's own mas- 
sive cl uster u Cen; both clusters have been suggested as possible remnant dwarf galaxy 
cores (IMeylan et al.l 120011). There has been some discussion about whether G001 contains 



a central black hole (IGebhardt et al.l 120021 ; iBaumgardt et al.l 120031 ; iGebhardt et al.l 120051 ) 



an d also some widely d iffering measurements of its structural parameters (see discussion 



m 



Barmbv et all 120021 ) . IGebhardt et al.l (l2005h used new HST/ACS observations of G001 



(obtained as part of GO-9767) to produce a surface brightness profile as part of an effort 
to constrain the mass of any putative centr al black hole, bu t t hose author s did not give 
values for the cluster's structural parameters. iMa et al.l (120071 ) fit iKingl (I1966T) models to the 



ACS d ata. Because of the importance of G001, we decided to re-analyze the IGebhardt et al. 
(120051 ) observations, which consist of six exposures totaling 41 minutes of integration time 
with the ACS HRC in the F555W filter. 

The HST images offered the possibility of detecting new clusters or confirming addi- 
tional cluster candidates other than the targets. All the M31 images were searched for 
additional clusters and inspected at the positions of cataloged cluster candidates. A number 
of additional cluster candidates appeared in the ACS images: B041, B056D, B061, B081D, 
B088D, B090, B102, B14 7, B162 and M027. We also (independently) re-discovered the new 
cluster B515 reported by Galleti et al. ( 20061 ). All of these candidates appear to be clusters 
except for B102, which is apparently a star@ Interestingly, B147, which had been cla ssi- 
fied as a star based on velocity dispersion measurements by iDubath fc Grillmairl (119971 ). is 
clearly a resolved cluster (see Figur e CD). This object is also associated with an X-ray source 
( ITrudolyubov fc Priedhorskyj 120041 ) . Figure [1] shows images of B147 and the highly elliptical 
cluster B088 demonstrating the power of HST imaging for resolving M31 globular clusters. 

To generate the largest possible sample of M31 globulars for fundamental plane analy- 
sis, we also include most of the clusters with structural parameters measure d from images 



taken with HST's Wide Fi eld and Planetary Cam era 2 (WPFC2), reported in lBarmby et al. 



( 120021 ). The clusters from iBarmby et al.l (120021 ) included in the present sample include 59 
objects with HST/WFPC2 photometry in the l/-band and good measurements of r and 



1 Galleti et al.l (|2006l ) independently came to the same conclusions about the nature of B056D, B102, and 
B162. 
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c. The central surface brightn ess values for thes e cluste rs have had an important error 
corrected. As a lso discussed in IMcLaughlin et al.l (l2007af) . the so ftware used to fit struc- 



tural models in Barmby et al. ( 20021 ) and also Harris et al. (bo02h did not return a model 
central surface brightness value; instead both papers reported fiyfl as the surface bright- 
ne ss of the centr a l pixe l of the PSF-convolved model image. As such, the values rep orted 



m 



Barmbv et al.l (120021 ) are syste matically too fa i nt. W e re-computed /iy,o from the [King 



(119661 ) c an d r n values reported in IBarmbv et al.l (120021 ) and the integrated V magnitudes 



reported in IBarmbv et al.l (120001 ) or IBarmbv fc Huchral (120011 ). The median surface bright- 
ness change is Ap,y,Q — —0.33 mag, with a large scatter. The maximum change is about 
1.5 mag, and there were several clusters, both bright and faint, whose central surface bright- 
nesse s were nearly unchanged in t he re-calculation. U sing the model-fitting software of 



both lMcLaughlin et al.l (j2007af ) and lBarmby et al.l (120021 ) on a small number of M31 clusters 
showed that both recovered the the same r and c values to within the uncertainties, so the 
two sets of measurements should be consistent. 

To summarize, the new observations include 15 clusters observed using the ACS in ei- 
ther the Wide Field or High-Resolution Channel, in both F606W and F814W filters; and 
19 clusters observed in the unfiltered mode using STIS. The observations of cluster G001 
with ACS/HRC in the F555W filter are not new but are re-analyzed here. There are there- 
fore 50 surface-brightness profiles for 34 distinct objects (one of the STIS clusters, B023, is 
also i n the ACS sample a nd has thus been observed three times). Adding the 59 objects 
from IBarmbv et al.l (120021 ). the sample of M31 clusters totals 93 objects, with some later 
omitted from the fundamental plane analyses (see § I3.3p. Estimates of the total size of 
the M31 GCS vary: one recent value is N gc = 460 ± 70 (IBarmbv fc Huchral l200ll ). Thus 



the HST sample comprises only a fraction of the total M31 population, but it is com- 
parable in s i ze to the full sample of non-core-collapsed Milky Way clusters considered by 
McLaugh lin (2 0001). and about 10% larg er than the sample of 85 Milky Way clusters analyzed 
by IMcLaughlin fc van der Marel ( 2005 ). We have not attempted to generate a complete 
sample of clusters observed with HST by searching the Ar chive for clusters serendipitously 
observed in data taken after the IBarmbv fc Huchral (120011 ) sample was compiled; however 
we expect that the present sample should be representative of M31 GCs particularly at the 
high-luminosity end. The exception is low-surface-brightness clusters for which our present 
analysis method is unsuitable: these include both faint clusters, some discussed in the fol- 
lowing section, and also the 'extended luminous' clusters discovered by iHuxor et al.l (120051 ) 
(see also §6.31) . 
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3. Data analysis 



The data analysis procedures used for these data are described in detail by lMcLaughlin et al 



( j2007al ). Here the procedures are briefly summarized, and some details specific to the M31 



observations noted. 



3.1. Surface Brightness Profiles 

The STScI-pipeline output drizzled images were used to measure cluster shapes with 
ELLIPSE in IRAF. Cluster center positions were fixed at values derived by centroiding and 
elliptical isophotes were fit to the data, with no sigma-clipping. A first pass of ELLIPSE was 
run in the usual way, with ellipticity and position angle allowed to vary with the isophote 
semi-major axis. Luminosity-weighted averages of ellipticity and position angle were deter- 
mined from the results. In the second pass of ELLIPSE, surface brightness profiles on fixed, 
zero-ellipticity isophotes were measured (this was required because we chose to fit circular 
models for both the intrinsic cluster structure and the PSF). For several low-density clusters 
(B290 and B423 observed with STIS; B081D, B088D, B515, and M027 observed with ACS), 
a sensible, monotonically-decreasing profile could not be obtained: no models were fit to 
these clusters. Table [1] reports the average ellipticities, position angles, and aperture mag- 
nitudes measured from the ACS and STIS images, and tabulates some additional reference 
data (metallicity, reddening, galactocentric distance) for all of the M31 clusters. For 11 of 
the clusters, there are no metallicity estimates in the literature, so we have simply assigned 
[Fe/H] = —1.2 ± 0.6 to all of the m (the mean and sta ndard deviation of the metallicity 



distribution of the M31 GC system; iBarmby et al.l 120001 ) . 



The raw output from ELLIPSE is in units of cts/s per pixel, which is converted to cts/s 
per square arcsecond by multiplying by 400 = (1 px/0'/05) 2 for WFC, 1600 = (1 px/0'/025) 2 
for HRC, or 387.5 = (1 px/0'/0508) 2 for STIS. To transform the ACS counts to surface 
brightness calibrated on the VEGAMAG system (from the ACS Handbook) 

/x/mag arcsec -2 = z — 2.5 log(cts/s/D") (1) 

where the zeropoints are z = (25.255, 26.398, 25.893, 25.501, 24.849) for HRC/F555W, WFC/F606W, 
HRC/F606W, WFC/F814W, and HRC/F 814W respect i vely. The STIS data were taken in 



unfiltered ( "clear" or CL) mode, for which lHarris et al.l (120021 ) derive 



Hv /mag arcsec" 2 = (26.29 ± 0.05) - 2.5 k>g(cts/s/D") (2) 



in the standard V bandpass. 
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We allow for the occasional over-subtraction of sky in the automatic reduction pipeline, 
and hence "negative" counts in some pixels, by working in terms of linear intensity in- 
stead of surface brightness in magnitudes. Given_| M@ jF6mw = 4.64, M Q ^ FS14W = 4.14, and 
M Q y/ F555W = 4.83, the calibration above yields 



If555\v/ Lq pc~ 2 




2.8766 x (cts/s/D") 


[HRC] 


lF606w/ Lq PC~ 


f-^ 


0.8427 x (cts/s/D") 


[WFC] 


lF606w/ Lq pC~ 2 


f-^ 


1.3418 x (cts/s/D") 


[HRC] 


LpsUW/ Lq PC~ 


f-^ 


1.2147 x (cts/s/D") 


[WFC] 


IfsUw/ Lq pc~ 2 




2.2145 x (cts/s/D") 


[HRC] 


Iv/Lq pc" 2 




1.1089 x (cts/s/D") 


[STIS] 



(3) 



Converting from luminosity density in L & pc 2 to surface brightness is done according 
to /isi4/mag arcsec -2 = 25.712 — 2.5 log(Isu/ L Q pc~ 2 ) for F814W; /i606/ m ag arcsec -2 = 
26.212-2.5 log(I 6O6 /L pc -2 ) for F606W; /^555/mag arcsec" 2 = 26.402-2.5 log(iy i5 55/£© pc~ 2 ) 
for standard V and F555W. 

The ELLIPSE profiles extended to R = 12 / /8 ~ 48 pc from all images of the sample 
clusters, and thus in the model fitting a constant-background term is included to allow 
for any errors in sky subtraction in the automatic pipeline. Each profile was inspected, 
and one or two very discrepant points in each of a handful of clusters were excluded by 
hand. Inspection of the data quality images for the WFC images showed that, due to an 
unfortunate choice of CCD gain, the cores of several bright clusters (B042, B063, B082, B147, 
and B151) were saturated, with raw intensities above a nominal limit of 70 cts s" 1 pix -1 
(28,000 cts s -1 arcsec -2 , or Ifgogw —23, 600 L & pc~ 2 ) all the way out to R ~ 2" = 7.5 pc (or 
~40 pix) in most cases. This is well beyond the FWHM of the WFC PSF in either F606W 
or F814W; see below. Intensities as high as this were not included in the model fitting, and 
thus the core parameters of these five clusters are constrained only very indirectly. None of 
the HRC or STIS images is saturated. 

Table [2] gives the final, calibrated intensity profiles for all clusters. These are not cor- 
rected for extinction, which is discussed below. The reported F606W- and F814W-band 
intensities are calibrated on the VEGAMAG scale, while the STIS data are on the standard 
V system. The final column gives a flag for every point, which can take one of four values. 
"BAD" indicates that the intensity value is deemed dubious because it strongly deviates 
from its neighbors or is obviously affected by nearby bright stars or image artifacts. "SAT" 
indicates that the isophotal intensity is above the imposed saturation limit of 70 cts s" 1 pix -1 



3 http : //www . ucolick . org/~cnaw/sun . html 
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"DEP" indicates that the radius is inside R < 2 pix = Cf!l and the isophotal intensity is 
dependent on its neighbors. The ELLIPSE output includes brightnesses for 15 radii inside 
2 pixels, but these are all measured from the same 13 central pixels and are clearly not sta- 
tistically independent. To avoid excessive weighting of the central regions of the cluster in 
the fits, only intensities at radii R m i n , R m i n + (0.5, 1.0, 2.0px), R > 2.5 px, were used. "OK" 
indicates that none of the above apply and the point is used in model fitting. 



3.2. Point-Spread Functions 



M31 clusters are clearly resolved with HST, but their observed core structures are still 
affected by the PSF. We chose not to deconvolve the data, instead convolving the structural 
models with a simple analytic description of the PSF before fitting. 

To estimate the PSF for the WFC, ELLIPSE was used (again with circular symmetry 
enforced) to produce intensity profiles out to radii of about 2" (40 pixels) for a number of 
isolated stars on a number of images, and combined them to produce a single, average PSF. 
This was done separately for the F606W and F814W filters. We originally tried to fit these 
with simple Moffat profiles (with backgrounds added), but found that a better description 
was given by a function of the form 



J PSF //o = [1 + (R/r ) 3 ] 



-P/3 



(4) 



The derived r and (3 values, and the implied FWHMs, are liste d in Table El Although the 
WFC PSF is known to vary over the instrument's field of view, iMcLaughlin et al.l (j2007al ) 
show that the model-fitting results are insensitive to this effect. 

For HRC data, the same functional form (eqB]) wa s fitted to the PSF images in F555W, 
F606W and F814W given by lAnderson fc King (120051 ) with the results also list ed in Table 
|3l In the case of STIS, a similar procedure was followed (see lHarris et al.ll2002l ) and it was 
found that a standard Moffat function, Jpsf oc [1 + (-R/Yo) 2 ] - ^ 2 , gives an adequate fit to 
the PSF, with parameters also listed in Table [3j These are the functions convolved with the 
structural models before fitting to the observed intensity profiles. 



3.3. Extinction and Color Corrections 



The effect ive wavelengths of t he ACS F606 W and F814W filters are A cfr = 5918 A and 



A e ff = 8060 A ( Sirianni et al.ll2005l ). so that from lCardelli et al.l ( 
and t4 814 ~ 1.8 x E(B — V) — both of which are also found by 



989 ) A mH ~ 2 . 8x E( B-V) 



Sirianni et all (120051 ) . This 
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implies, rather conveniently, that E(F606W — F81AW) ~ E(B — V). Individual reddening 
estimates — the sum of foreground and M31-internal reddening — are available for most of the 



clusters in the sample, d erived using the procedure described by iBarmby et al.l (120001 ) . The 
Burstein fc Heilesl ( 11984 ) foreground value of E(B — V) = 0.08 mag is adopted for clusters 
without individual measurements. 



Sirianni et al.l (120051 ) give transformations from ACS F606W and F814W magnitudes 
to standard V and I for both WFC and HRC measurements, including linear and quadratic 
dependences on de-reddened (V — J) color (their Tables 22 and 23). However, measured 
(V — I) is not available for 5 of the 15 clusters in the present ACS sample, and thus we 
chose instead t o measure ape rture (F606W — F81AW) colors from the current data, and 
manipulate the ISirianni et al.l relations to find (V — FQ0QW) and (V — F81AW) indices. 



Comparing the various "observed" and "s ynthetic" transforma tions from intrinsic FQ06W 
magnitudes and (V — I)q colors to V mags in lSirianni et al.l (120051 ). it quickly becomes clear 
that there are systematic differences between them, at the level of ~0.07 mag in the color 
range 0.8 < (V — 1) < 1.2 appropriate for old clusters with sub-solar metallicities. However, 
a simple linear relation, 



{V - F6O6WO0 = -0.038 + 0.236(V - I\ 



(5) 



both gives a good estimate of the average of the ISirianni et al.l relations (for VEGAMAG 
F606W magnitudes ) and agrees very well w ith the predictions for old clusters in the population- 
synthesis models of lMarastonl Jl998l . l2005h . [VEGAMAG (V - F606W) and (V - F8UW) 
indices were kindly computed for us, as functions of cluster age and metallicity, by C. Maras- 
ton.] 

To go further, VEGAM AG F8UW magnitud es must be related to /. All of the WFC 
and HRC transformations in ISirianni et al.l (120051 ) imply that (F81 4W — I)q ~ to within 
±0.01 mag or so, while the population-synthesis models of lMarastonl have a roughly constant 
(F8UW-I) ~ 0.04 for 0.8 < (V- 
© yields 

(V - F6O6WO0 



I)o < 1.2. Substituting (F8UW-0.0A) for / in equation 



and then trivially, 



-0.04 + 0.31(F606W/ - F8UW\ 



{V - F814WOo = -°- 04 + 1.31(F6(W - F8UW) Q 



(6) 



(7) 



both with estimated rms scatter of about ±0.05 mag. These conversions were applied 
to both WFC and HRC data, after correcting the calibrated F606W and ^8141^ mag- 
nitudes/intensities for extinction. 
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Calibrating the other observations onto the standard V scale was more straightforward. 
The ACS /HRC F555W filter is similar to the standard V filter, and the data for G00 1 were 
calibrated using the extinction ratio and transformation given by lSirianni et al.l (120051 ) (their 
Tables 15 and 23). The dat a for the 19 STIS clusters are al ready calibrated on t he sta ndard 
V scale lHarris et all (boOSZh . as are the WFPC2 data from barmbv k Huchral (boOlh . The 
reddening values used to correct for extinction are recorded in Table [TJ 

Given the color transformations above, the measured ACS aperture colors plus the 
assumed reddening values allow us to predict a reddened V — I color to be compared with 
ground-based aperture colors. For the 10 of 1 5 ACS-observed clus ters with (V — 7) colors 
(ground-based, mainly from the compilation in lBarmby et al.ll2000l ). the agreement between 
predicted and observed (V — I) colors is very good: the mean offset is 0.02 ± 0.02 mag, with 
a median absolute deviation of 0.06 mag. The cluster B023 has been observed in all three 
ba ndpas s es dis cussed here: F606W and F814W with ACS/HRC, and V with STIS. The fits 
of iKingl (119661 ) models to these three independent profiles in §|5] yield extinction-corrected 
total magnitudes of F8UW = 12.09, FQ06W = 12.88, and V = 12.99, with uncertainties 
of about ±0.02 mag in each. The implied global (F606W — F81AW)o = 0.79 is essentially 
identical to the aperture color computed from the measurements in Table [Tj and the fitted 
global (V — F606W) and (V — F814W)q are within 0.1 mag of the predicted values- 
suggesting a level of agreement similar to that between (V — I) prc d and (V — I) bs more 
generally. 

The cluster B041 is somewhat unusual. Its integrated colors suggest that it is perfectly 
average, compared to the rest of the M31 cluster system. However, the ACS image shows 
that a very bright, red star about 0'/5 from the cluster center is responsible for the redness 
of its (F606W — F814W) and (V — I) aperture colors. Both the innermost intensities 
from the ELLIPSE profiles and the total cluster F606W and F814W magnitudes from 
model fitting (in which the star is "smoothed over" in some sense), imply a de-reddened 
(F606W — F81AW)o — 0.4 — much bluer than any other cluster in the sample. In fact, 
the population-synthesis model predictions of C. Maraston suggest that the oldest possible 
age for such a blue cluster is < 1.5 Gyr. Because of its red integra ted color, this cluster is 
not included in any of the 'massive young M31 cluster' samples of iFusi Pecci et al.l (120051 ) 
- showing that integrated properties are not always reliable guides to the true nature of 
extragalactic globular clusters. There are eight clusters in the W FPC2 sample which ar e 
members of 'blue luminous compact cluster' samples 'A' or 'B' in IFusi Pecci et al.l ( 120051 ). 
These are suspected of being young clusters by virtue of very blue colors [(B — V) < 0.45] or 
strong H/3 absorption. These clusters and B041 are not included in any of the fundamental 
plane correlations in §6.3j a total of 84 M31 clusters are considered there. 



- 11 - 



4. Models 

Converting luminosity-based parameters from the model fits into mass-based quantities 
requires knowledge of mass-to-light ratios. Velocity-dispersion measurements do exist for a 
fraction of t he sample fsee ^5. 3D but t he measurements are far from complete or uniform. 
Instead, as in lMcLaughlin et al. (b()07ah . we use po pulation-synthesis models t o predict mass- 
to-light ratios for the full sample, u s ing th e code of iBruzual fc Charlotl (120031 ) and assuming 
the (disk) stellar IMF of IChabrierl (120031 ) . These can then be used to produce (for exam- 
ple) predicted velocity dispersions that can be compared to current and future spectroscopic 
data. The values of Y pop for an assumed age of 13 Gyr are discussed further and tabu- 
lated in §5.21 and shown as a histogram in Figure [2J For comparison, also shown is the 
(re-normalized to the same sample size) d istribution of mass-to-light ra t ios fo r 85 Galac- 
tic globulars with parameters cataloged by iMcLaughlin fc van der Marell (120051 ). calculated 
using the same population-synthesis model code with the same assumed stellar IMF and 
for the same assumed age of 13 Gyr as in the present analysis. The distribution peaks at 
a slightly higher value of Ty° p in M31 vs. the Milky Way, which simply reflects a slightly 
higher average GC metallicity in the sample of M31 clusters. 



The cluster structural models used in the fits are described in d etail bylMcLaughlin et al. 



( 12007al ). The three models considered here are the often-used King! (119661 ) single-mass, 
isotropic, modified isothermal sphere; an alternate modificatio n of a s i ngle-m ass, isotropic 
isothermal sphere, based on a model originally i ntroduc e d by IWilsonl (119751 ) for elliptical 
galaxies; and t he Sersid ( 19681 ) or R l / n model. IWilsonl ( 19751 ) models are similar to the 
standard iKingi (119661 ) models, in that their cores are nearly isothermal, but their envelope 
structures are relatively more extended — a feature that has been found to better describe 
the density profile s of old globulars and young massiv e clu sters in the Local Group and 
in NGC 5128. See IMcLaughlin k van der Marell (bood ) and IMcLaughlin et all J2007ah for 
more detailed discussion and for general comparisons bet ween the two types of models. Two 
additional models , 'power -law' profiles and the analytic IKingi (119621 ) model, are found by 
McLaughlin et al.l (j2007al ) to provide little or no information beyond that contained in the 
first three models, and are not discussed here. 

Fitting the structural models to the data involves first selecting a value for tq, the scale 
radius@ and computing the dimensionless model profile J mo d = I mo d/Io- The model is then 
convolved with the appropriate PSF having parameters taken from Table [3] to yield the 
product I^ od (R\ro). The fitting procedure allows for a varying scale radius r and non-zero 



4 The definition of r$ varies between models: for Sersic Jl968) models, I(rp) _ 0.5 Jn, while for Kine ( 1966f) 
and lWilsonl (| 19751) this relation is only approximate. See IMcLaughlin et al.1 (|2007a[) for details. 
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background Jbkg, and so minimizes 

2 

(8) 



* 2 = £ 



^obs(-Rj) - Iq X -fmod(-^il r o) - 4kg 



a? 



for the measured intensity profile I b s (Ri) and uncertainties Oi of each cluster in Table El 

Uncertainties on the fitted and derived model parameters are estimated in the usual way, 
from the range of their values in all models for which x 2 is within some specified distance of 
the absolute minimum (e.g., A% 2 < 1 for 68% confidence intervals). Because the formal error 
bars estimated by ELLIPSE for the isophotal intensities are artificially small, the best-fit 
Xmm can be exceedingly high (3> iVpts) even when a model fit is actually very good. This 
would result in unrealistically small estimates of parameter uncertainties. We chose to re- 
scale the x 2 f° r an fitted models by a common factor chosen to make the global minimum 
Xmin = (-Wpts — 4), where iVp ts is the number of points used in the model fitting. Under 
this re-scaling, the global minimum reduced x 2 P er degree of freedom is exactly one. This 
prescription for deciding fit uncertainties does not affect the identification of the best-fitting 
model itself, since the relative sizes of the uncertainties on individual data points are not 
changed. 



5. Results 



5.1. Model-fitting 



Figure [3] shows example surface brightness profiles and model fits for four M31 clusters 
chosen to span the observed luminosity range. The observed data plotted have had the fit- 
ted Jbkg subtracted and then been converted to (extinction-corrected) surface brightnesses. 
Points with intensities below the subtracted background in the upper panel are represented 
here as solid points placed on the lower x-axis, with error bars extending upwards. The Fig- 
ure shows that the M31 GCs are well-resolved by the HST imaging, and also demonstrates 
the core saturation for the brightest clusters exemplified by B151. For these particular clus- 
ters there is r elatively little difference in the goodness-of-fit of the three different models. As 
expected, the IWilsonl (119751 ) models, with t heir extended halos, favor somewhat lower back- 
ground levels, and the 'cuspier' ISersid (119681 ) models favor higher central surface brightnesses. 
The clusters have projected half-light radii Rh = 2 — 4 pc an d King-model conce ntrations 
in the range 1-1.5, quite typical values for GCs in both M31 (IBarmby et al.ll2002l ) and the 
Milky Way dHarrislll996h . 



Table H] summarizes the basic ingredients of all model fits to the full sample. The first 
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column in this table gives the cluster name; the second, the detector/filter combination from 
which the density profile was derived; the third, the extinction in the F606W band; the 
fourth, the color term to transform photometry from the native bandpass of the data to the 
standard V scale (see §3.3p . and the fifth column lists the number of points in the intensity 
profile that are flagged as OK in Table [2] above, and thus were used to constrain our model 
fits. Subsequent columns in Table H] cover three lines for each cluster, one line for each type 
of model fit. Column (6) identifies the model. Column (7) gives the minimum unreduced x 2 
obtained for that class of model (without the re-scaling applied for uncertainty estimation). 



King 


1966 


and 


Wilson 



Column (10) gives the concentration c = log(r t /Yo), or the ISersid index n of the best fit. 
Column (11) gives the calibrated and extinction-corrected central surface brightness in the 
native bandpass of the data. Column (12) gives the logarithm of the best-fit scale radius r 
in arcsec (see §1]), and Column (13) gives the logarithm of r in units of pc (obtained from 
the angular scale assuming D = 780 kpc for M31). E rror bars on all these parameters were 



defined in the same way as iMcLaughlin et al.l (j2007al ). 



5.2. Derived quantities 

Table contains a number of other structural cluster p roperties derived fr o m the basic 
fit parameters. The details of their calculation are given in IMcLaughlin et al.l (j2007al ); the 
contents are: 



1. log r t = c + log r , the model tidal radius in pc (always infinite for ISersid models) 

2. log R c , the projected core radius of the model fitting a cluster, in units of pc 

3. log Rh, the projected half-light, or effective, radius of a model: that radius containing 
half the total luminosity in projection, in units of pcj^l 

4. log (Rh/R c ), a measure of cluster concentration 

5. log Iq , the logarithm of the best-fit central (R = 0) luminosity surface density in the 
V band, in units of L & y pc~ 2 

6. log 7n, the logarithmic central (r = 0) luminosity volume density in the V band (for 
Sersid models the density at the 3-D radius r ) , in units of L Q pc~ 3 



5 Calculations show that the three-dimensional half-light radius = 4/3Rh- 
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7. log Ly, the logarithm of the total integrated model luminosity in the V band 

8- V t ot = 4.83 — 2.5 log (L V /L Q ) +5 log (D/lOpc) is the total, extinction- corrected appar- 
ent V-band magnitude of a model cluster 

9. log Ih = \og(L v /27iRl) is the V"-band luminosity surface density averaged over the 
half-light /effective radius, in units of L Q y pc~ 2 . 



The uncertainties on all of these derived parameters have been estimated (separately for 
each given model family) following the y 2 re-scaling proced ure described ab ove. The re- 
computed values for the 59 clusters from lBarmby et al.l (120021 ) . for [King (119661 ) models only, 
are reported at the end of this Table. 

Table [6] next lists a number of "dynamical" cluster properties derived from the structural 
parameters already given plus a mass-to-light ratio. The first two columns of this table 
contain the name of each cluster and the combination of detector /filter for the observations. 
Column (3) lists the mass-to-light ratio, in solar units, adopted for each object from the 
analysis in §H The values of Ty° p assume an age of 13 Gyr for all clusters. The error bars 
on T y op in Table [6] include a ±2-Gyr uncertainty in age, as well as the previously tabulated 
uncertainties in [Fe/H]. The remaining entries in Table [6] are, for the best fit of every model 
type to every cluster: 



1. log M tot = log T 



pop 



log Ly, the integrated model mass in solar units, with log L\ 



taken from Column (10) of Table [5] 

2. log Eb, the integrated binding energy in ergs, defined through Ej, = —(1/2) J Q r * 47rr 2 p0 dr 

3. log S , the central surface mass density in units of M pc~ 2 



4. log p , the central volume density in units of M pc 3 (for ISersid models the density 
at the 3-D radius r ) 

5. log S/j, the surface mass density averaged over the effective radius Rh, in units of 

M pc~ 2 

6. log cr p n, t he predicted line-of-sight velocity dispersion at the cluster center (or at R c 
for ISersid models), in km s -1 



7. log f eS c,o, the predicted central "escape" velocity in km s 1 

8. log t r h, the two-body relaxation time at the model projected half-mass radius, in years 
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9. log/o = log [po/(27i(Jc) 3 ^ 2 ] , a measure of the model's central (at r for ISersid models) 
phase-space density or relaxation time, in unit s of M M pc~ 3 (km s -1 )~ 3 . For f n in t hese 



units, and t rc in years, logt rc ~ 8.28 — log/o ( McLaughlin Sz van der Marel 20051 ) 



The uncertainties in these derived dynamical quantities are estimated from their variations 
around the minimum of x 2 on the model grids, as above, combined in quadrature with the 
uncertainties in the population-synthesis model T 



pop 



Finally, Table [7| provides the last few parameters required to construct the fundamental 
plane of globular clusters in M31 under any of the equivalent formulations of it in the 
literature (see §6.31 below). The last t hree columns of Table [7] are modified versions of the 

parameters of iBender et al.l (119921 ). who found that early- type galaxies and (separately) 
globular clusters defined the edges of planes in k space. We define the parameters using 
mass rather than luminosity surface density, as this is more u seful for comparing globulars 
to younger clusters and galaxies (see McLaughlin et al.ll2007bl ). and emphasize this by using 
the notation K m ^. 

(log (7*0 + log R h )/y/2 

(log oj )0 + 2 log Z h - log R h )/VE (9) 
( lo g °"p,o - lo £ - log Rh)/V3 

In calculating « m ,i, «m,2> an d K m,3 for Tabled the cr Pi o predicted in Table [6] (Column 7) 
by adoption of population-synthesis mass-to-light ratios has been used, and Y* h evaluated 
by adding log Ty° p to log Ih from Column (12) of Table [5l (This means is therefore 
independent of the assumed mass-to-light ratio.) Rh is also taken from Table [5] but put 
in units of kpc rat her than pc, for compatibility with the galaxy-oriented definitions of 
Bender et all Jl992h . 



5.3. Predicted Velocity Dispersions 

One use for the material presented in the previous subsections is to predict observable, 
line-of-sight velocity dispersions averaged over circular apertures of any radius, which can 
then be compared directly to extant and forthcoming spectroscopic data on the clusters in the 
sample. Ultimately such comparisons will be useful not only for independent assessments of 
the self-consistency of the structural model fits, but also to evaluate how close the dynamical 
mass-to-light ratios in these old globulars come to the expectations of population-synthesis 
models with no dark matter. 

Given any of the models with a fitted r and c or n for any cluster, solving Poisson's 
and Jeans' equations and projecting along the line of sight yields a dimensionless velocity- 
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dispersion profile, a p = a p (R)/<7 P)0 as a function of projected clustercentric radius R. The 
density-weighted average of within a circular region of radius i? ap then gives a dimension- 
less aperture dispersion for the cluster. Finally, normalizing with the predicted cx Pi o from the 
fits in Table El we can predict an observable cr ap (i? ap ) for every model for every cluster, and 
for any specified _R ap . Table [8] gives the results of these calculations for each of 5 apertures: 
i? ap = Rh/8, Rh/4:, Rh,4: Rh, and R ap = 8Rh- For a distance of 780 kpc to M31, a typical 
GC Rh = 3-4 pc translates to roughly O'.'8-l'.'l, so the range of apertures over which cr ap 
is calculated should be a good match to reasonable observations. The velocity dispersion 
inside other apertures can be obtained by interpolation. 

These predicted velocity dispersions can already be compared in a limited way to ob- 
servations. Measured velocity dispersions for 2 9 M31 GCs, o f which 22 are included in th e 
current HST sample, a re av ailable from fro m I Cohen! (120061 ). iDubath fc Grillmairl (119971 ). 
Djorgovski et al.l (119971 ) , and iPetersonl (119891 ) . These are summarized in Table [9j Figure H] 
compares observed and model-predicted velocity dispersions, where the model predictions 
are derived from interpolation on the values in Table [8] to the aperture size reported for 
the observations. Aperture sizes were computed as R ap = yjlw/ir where I and w are the 
reported spectrograph slit dimensions. 

Figure H] shows that, except for a constant ratio, there is generally good agreement 
between observed and predicted velocity dispersions (the exception, cluster B037, is discussed 
further in §5.51) . This justifies our use of population synthesis model mass-to-light ratios 
to derive masses and other properties such as binding energy. This procedure gives us a 
much larger sample of clusters to work with for a wider variety of analyses than would 
be possible using only direct velocity dispersion measurements. The median ratio between 
observed and predicted velocity dispersions is o" bs/o" pre d = 0.85, with an inter-quartile range 
±0.15. This corresponds to an ratio between dynamical and population-synthesis-derived 
mass-to-light ratios of x^ yn /Ty° p = 0.73 ± 0.25 , consistent with the value for this ratio of 



0.82 ± 0.07 found bv McLaughlin fc van der Marell (120051 ) for Milky Way and old Magellanic 
Cloud clusters. The lower dynamical mass-to-light ratios compared to those derived from 
population-synthesis models is consistent with with theoretical expectations that old Milky 
Way and M31 GCs have lost low-mass stars through evaporation, an effect not included in 
the models. 



5.4. Comparison of Fits in F606W and F814W 



Comparing model fits to the same cluster observed in different filters allows assessment 
of the systematic errors and color dependencies in the fits. Figure compares the parameters 
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derived from fits to the 15 ACS clusters observed in both F606 W and F814W . In g eneral 
the ag reement is quite good, with somewh at higher scatter for the IWilsonl (119751 ) and ISersic 
( 119681 ) model fits than for the iKingi (119661 ) fits. There are a few unsuccessful fits with very 
high c oncentrations, small sizes and extreme l y brig ht central surface brightnesses: the iKing 
jl966h model fit to B082 in F606 W; IwilsorJ Jl975h model fit to B042 in F606W and B063 
in F814W; and the ISersid Jl968h model fit to B147 in F814W. These are all clusters in 
which the central core data are not used due to saturation; fortunately the data in the 
other filter produce a more reasonable-looking fit. For all three models, the central V-band 
surface brightnesses as fit to the F814W data are slightly fainter than those from the F606W 
data, while the total model luminosities are essentially the same. The latter implies that our 
estimated (V — F606W) and (V — F814W) colors in §3.3l are correct, while the former implies 
that the cluster centers are slightly bluer than the average cluster color. This could indicate 
the presence of blue stragglers or blue horizontal branch stars in these massive clusters, but it 
could also be due to the differing levels of saturation. A broader range of data will be needed 
to investigate the issue of color gradients in M31 GCs more thoroug hly (for a discussion of 
color gradients in Milky Way clusters, see iDjorgovski fc Piottolll993l ). 

Because the model fit results in the two ACS bands are quite similar, and also because for 
the STIS and WFPC2 data we have only l/-band data, we do not consider the F814W model 
fits in the following analyses. Parameters derived from the model fits given in Tables [5H8] 
refer to the F606W band observation s only. The exception is B082, where we substitute the 
F814W fit results for the[King| Jl966h fit. 



5.5. Comparison to previous results for individual clusters 



see 



The famous M31 c l uster G001 has discrepant published values of its half-light radius 



Barmbv et al. 



2002 



in this paper. iMeyl an et al 



a projected Rh = 2'.'8, while 
Rh = 0"70fj Our 
light radius R 



Ma et al.ll2006bl ). which was a major motivation for re-analyzing it 



(2001 ) reported a th ree-dimensional = 3'/ 7, which converts to 



Barmbv et all (120021 ) measured R h = 0"82 and lRich et all (11999 ) 



Kinej Jl966h model fit s to th e ACS/HRC F555W data resulted in a half- 
/ /85, and our IWilsonl (119751 ) fit to the same data yielded Rh = 1"2. The 



Meylan et al.l (120011 ) value is clearly the outl ier, likely because it is for the half- mass, rather 
than t he ha lf-light radius, and because the iMeylan et al.l (120011 ) analysis used multi-mass 
King! (119661 ) models. The surface brightness profile from the HRC observations is almost 



6 lLarseij (|200lf ) pointed out that the integrated photometry given by iMevlan et al.l (|200ll ) implies a 
half-light radius of about 1'.'18, much closer to other values. 
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servations including those used by lRich et al 



identical to that presented in lMevlan et al.1 (20011). which combined data from WFPC2 ob 



(119961 ) , so we conclude that the physi cal model 



assum ptions used are responsible for the differences between authors. Recently, iMa et al. 



( 120071 ) also analyzed the ACS/HRC data used here: t hese authors derived a higher concen- 
tration c = 2.01 ± 0.02 compared to our [King (119661 ) model value of 1.77, and a value of 
Rh = 1'.'73 ± 0'.'07 about twice as large as ours. This appears to be the result of weighting 
the outermost points in the surface brightness profil e more he avily, yielding a larger tidal 
radius and higher concentration. The dashed line in IMa et all 's Figure 3, which represents 
their model fit with the outermost surface brightness points excluded, is quite similar to our 
model fit, and we believe this to be more robust particularly in the better-defined inner part 
of the cluster. 

The heavily- reddened M31 cluster B037 has also been the subject of intensive s tudy- 



Ma et all (j2006bl ) analyzed the same F606W image of B037 discussed here, fitting a [King 
( 19621 ) model to a surface brightness profile made from a PSF-deconvolved image. They de- 
rived a core radius ro = 0'.'72, half-light radius = l'.'ll, concentration c = 0.91, and central 
surface brightness /i(0) = 17.21 mag arcsec -2 in F606W [using our values for extinction and 
V — F606W / color, this becomes \xy$ = 13.49.] Our model fits generally result in a some- 
what highe r concentration and smaller core radius [r = 0'.'56, = 1'.'09, c = 1.23 for the 



King! (119661 ) model]. This relatively small change would not likely have affected IMa et al 
conclusions about the nature of this cluster. 



Ma et all (j2006al ) predict a velocity dispersion of 72 km s _1 for B037, on the basis of a 
mass estimate of ~ 3 x 1O 7 M . However, these M and a values are based on a population- 
synthesis M/Ly ratio applied to an intrinsic cluster luminosity obtained by assuming a 
distance of ~ 950 kpc to M31, instead of our adopted 780 kpc (and a large E(B — V) 
close to the one we have used). Thus, comparison of their numbers wit h ours req uires 



first multiplying their derived mass by (780/950) 2 = 2/3. In addition, the IMa et al 



mass 



and velocity dispersion are actually calculated for what the cluster would have been at an 
age of 10 Myr. Howeve r , the same population-synthesis models that they (like us) used — 



i.e., 



Bruzual &: Charlotl (120031 ) — show that the cumulative effects of supernova explosions, 



massive-star winds, and AGB mass loss over 12-13 Gyr lead to a reduction of the 10- 
Myr total cluster mass by ano ther factor of ~ 2 at its current, old age (see also, e.g., 
Fioc fc Rocca-Volmerangelll997l ). Therefore, the present-day mass of B037 implied by the 
work of IMa et al is in fact ~ 1 x 1O 7 M , in good agreement with our results here (Table [6]). 



The high velocity dispersion predicted for B037 by IMa et all is a 'global' (not central) 
value following from the virial theorem assuming both their high mass and a rough estimate 
of the cluster half-light radius (~ 2.6 pc in their analysis), which is significantly smaller than 
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that in our more accurate profile fitting (Rh = 3.9 pc; e.g., Tabled column 5). Scaling to a 3- 
times smaller present-day mass, as above, and to a ~ 1.5-times larger Rh, the expected global 
velocity dispersion following self-consistently from Ma et al. is 72 km s" 1 x a/ (l/3)/1.5 = 
34 km s _1 , which is then entirely consistent with our predicted cr ap inside the largest apertures 
in Table E 



Cohenl (120061 ) has measured a velocity dispersion of 19.6 ± 3.5 km s -1 for B037 within 
an aperture of _R ap ~ 1.6 arcsec. Our predicted dispersion within this aperture follows 
from interpolation on Table El 38.3 ± 2.7 km s . This factor-of-two discrepancy (which 
implies a factor-of-four discrepancy between the true cluster mass and that estimated from 
our de-reddened total luminosity and 13-Gyr old population-synthesis M/L v ratio) could 
be resolved by lowering the mass-to-light ratio, but reducing Ty° p by a factor of 4 would 
require imply an age for B037 of only a few Gyr, in conflict with the SED-fitting results. 
Reducing the reddening of B037 to E(B — V) ~ 0.85 would yield a total l uminosity in better 
agreement with the measured a v , but this would also conflict with th e iMa et all (I2006ah 
SED-fitting. The only viable resolution appears to be the suggest ion by Co ienl (2006 ) that 



the dust lane projected across the face of the cluster (reported by IMa et al 



2006bJ) distorts 



measurements of extinction from integrated photometry. This should be testable with color 
profiles and star-count analysis from the ACS images. Such analysis is beyond the scope of 
the present paper, and the properties of B037 computed using E(B — V) = 1.36 are used in 
our subsequent analysis. 



5.6. Quality of Fit for Different Models 



McLaughlin fc van der Marell (120051 ) and iMcLaughlin et al.l (j2007al ) give detailed dis- 
cussions of the differences between various cluster models and their applicability to GCs in 
galaxies other than M31. Which structural model best represents the M31 clusters? Figure[6] 
compares the \ 2 values for the model fits as a function of cluster size Rh and total (model) 
luminosity Ly. The difference between models is parameterized as 



A = (x 



alternate 



Xk ins )/{x. 



alternate 



mg/ 



(10) 



which is independent of the \ 2 re-scali ng; significan tly better fits of one model over another 
are s ignaled b y [A > 0.2. In gener al, iKingl (119661 ) models fit the M31 cluster data better 
than IWilsonl (119751) or ISersid (119681 ) models. Cluster size does not strongly affect A, but 
luminosity does: King ( 1966h models ar e more strongly preferred for more luminous M31 



clusters. The preference for IKingl ( 119661 ) models for M31 clusters differs from results using 
the same profile-fitting software for for clusters in the Milky Way, Magellanic Clouds, and 
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2007af ). where 



Fornax dSph (IMcLaughlin fc van der Marelll2005l ) and NGC 5128 cluste rs (McLaughlin et al. 



Wilson! Jl975T models fit as well as or better than than iKingl (jl966l ) models. 



The results for M31 clusters do not appear to depend on the image saturation for some 
clusters — the clus ters observed with ACS/HRC, where saturation was not a problem, were 
also better-fit with iKind (119661 ) models. The fitting procedure returned higher background 
levels for M31 clusters than those in NGC 5128, but the background value does not correlate 
with A. 



Figure [7] compares the values of cluster parameters derived from lKingi (119661 ) and I Wilson 



( 119751 ) model fits. The parameters do not clearly vary with goodness-of-fit, and even the 
major outliers are not always much better fit with one model or the oth er. T he mean 



offsets between parameters derived for the same cluster from iKingl (119661 ) and IWilsonl (119751 ) 
models are comparable to the rms scatter, and only slightly larger than the typical parameter 
uncertainties: 5fi Vj0 = 0.04 ± 0.18, 8(\ogR h ) = 0.14 ± 0.13, 5{logL v ) = 0.08 ± 0.06, and 
5(logR r ) = .03 ±0.02. Small-scale systematic differences between the models (for example, 



Wilson! (119751 ) models, which have more extended halos, always return fainter central surface 
brightnesses /zy,o than IKingl ( 1966 ) models) do not affect the global cluster parameters, such 



King 



King 



as the binding energy a nd tota l lumin osity Ly . These are very similar in all three of 
(1966), IWilsonl (119751) and ISersid (119681 ) models. We suspect that the preference for 
( 119661 ) over IWilsonl (119751 ) models is due to some more subtle feature of the observational 
data that we have not yet been able to isolate, rather than an intrinsic difference between 
M31 clusters and those in oth er galaxies. H owever, since a majority of the M31 clusters have 
parameters derived only with IKingl (119661 ) model fits, we use these models as the basis for 
the fundamental plane analysis presented in §6.31 We do not expect the following discussion 
to depend on the model choice. 



Discussion 



We now combine the M31 GC parameters newly derive d in the previous section with 
those for clusters in the Milky Way and Mag ellanic Clouds (IMcLaughlin fc van der Marel 
20051 ) and NGC 5128 (IMcLaughlin et al.ll2007al ) to form a large sample of GCs that have been 
analyzed in a nearly homogeneous way. The sample comprises 291 clusters in six galaxies, 
an unprecedentedly large sample for defining the GC fundamental plane. In the following 
analysis, cluster luminosities are derived from the model fits to the surface photometry rather 
than integrated or aperture measurements. For the M31 and NGC 5128 clusters, model and 
measured luminosities are found to have excellent correspondence, while for Milky Way 
and Magellanic Cloud clusters, discrepancies between model and measured luminosities are 
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often attributable to the measurem ent aperture being smaller than the clusters' full size 
( [McLaughlin fc van der Marell 120051). Where M31 clusters are compared with NGC 5128 
clusters from McLaughlin et al.l (I2007af) , the 2 7 clusters omitted from that analysis are also 
omitted here. As discussed above, iKingl (119661 ) model parameters are adopted for clusters in 
all galaxies. 



6.1. Ellipticity distribution 



There are several possible explanations for the elongation of globular clusters: lLarsen 



(120011 ) lists elongation as possibly resulting from galaxy t ides, internal rotation, cl uster 
mergers (crossing times in star cluster complexes are short; iFellhauer fc Kroupal 120051 ) and 
'remnant elongation' (due to accretion?) from some clusters' former lives as dwarf galaxy 
nuclei. The latter two mechanisms would seem to be applicable to only a small fraction of 
clusters, probably the more massive one. They could be related to the observations that 
the most luminous young and old LMC clusters tend to be more flattened, and that the 
brightest globular clusters in both the Milky Way a nd M31 are also the most flattened 
(Ivan den Bergh fc Morbeylll984j ; Ivan den Bergdfl996[ )Fl Measurement of GC ellipticities is 
an uncertain business particularly for faint clusters: the shape of outer isophotes is affected by 
the galaxy background light and low signal-to- noise. For the extragalactic clusters measured 
with ELLIPSE, very low values of ellipticity are disfavored because the ELLIPSE algorithm 
diverges at e = 0. With these caveats in mind, it is still worthwhile to explore ellipticity 
distributions and correlations with other parameters for clues to the origin of GC flattening. 

Testing the idea that GC elongation is due to galaxy tides can be done to some extent 
by comparing the ellipticities of clusters in different galaxies, and thus different dynamical 
environments . Figure [8] shows the distribution o f ellipticities for cl usters in the Milky Way 
dHarrijll996h. M31 (this work), and NGC 5128 (Harris et al.ll2006h . Compared to Figure 7 
of lliarris et al. ( 2002 ). the number of M31 clusters is about 10% larger (although the sample 
is somewhat different and includes more very luminous clusters); the number of NGC 5128 
clusters is about 3 times larger and the distribution of their ellipticities somewhat more 
skewed to lower values. The distributions of ellipticities for M31 and NGC 5128 are not 
statistically different; both differ from the Milky Way distribution in having few very round 
clusters (likely an artifact of the measurement technique). The bottom panel of the Figure 
shows the distribution of ellipticity with galactocentric position: no correlation is evident. 



7 As discussed in §5.5[ there is some controversy about whether G001 or B037 is the most luminous M31 
GC, but both are quite elliptical. 
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From this comparison there is no clear evidence for major differences in ellipticity distribution 
between the three different galaxies' GCs, and thus no evidence that the overall galaxy 
environment is a major factor. 

If cluster ellipticities are caused by internal processes such as rotation or velocity 
anisotropy, then relaxation th rough dynamical evolution should act to reduce any initial 



flattening (IFiestas et al.ll2006l ). More dynamically-evolved clusters might be expected to be 
rounder (although the relaxation time is long in the outer regions of clusters which heavily 
influence measured ellipticity). Figure M shows ellipticity as a function of luminosity and 
half-mass relaxation time for clusters in the three galaxies. A mild systematic decrease in 
e with increased L is visible, though with considerable scatter; measurements of very large 
ellipticiti es for clusters w i th low luminosities are more likely to be spurious. While the obser- 



vation of Ivan den Berghl (119961 ) that the most luminous clusters in M31 and the Milky Way 
are highly elliptical remains true, it is clear that all three galaxies also have very luminous 
clusters which are quite round: did these bright, round clusters simply have less angular 
momentum at formation? Or are they the 'true' globular clusters while the more flattened 
luminous clusters are remnant dwarf galaxy nuclei somehow flattened while being removed 
from their original galaxies? Further modeling of the dwarf galaxy stripping process might 
help to answer this question. No correlation of e with the more important dynamical quan- 
tity t r h is evident: the data do not appear to favor the contention that more evolved clusters 
are rounder. While primordial angular momentum still seems to be the most straightforward 
explanation for GC ellipticity, the observed distribution appears to be due to a number of 
factors. 



6.2. Correlations with position and metallicity 



Structures of GCs in the Milk y Way have been shown to be largely independent of galac- 



tocentric distance and metallicity ( Djorgovski fc Meylan 1994 : Djorgovskjl 
200oh . except for the correlation of half-light radius with R ? r first noted by 



(1991J). Studies of cluster structures in nearby galaxies (IBarmby et al. 



1995 



McLaughlin 



van den Bergh et al 



2002 



Harris et al. 



20021 ) have found similar results, with the exception of a mild correlation between clus- 
ter half-light radius and metallicity: more metal-rich clusters tend to be smaller. The 
same e ffect is seen in more dis tant galaxies where the clusters' structures are barely re- 
solved. iBrodie fc Straderl (120061 ) discuss several possible explanations and conclude that the 
Rh — [Fe/H] correlation is most likely a consequence of the correlation between Rh and R gc 
and the tendency of metal-rich GCs to be located at smaller R gc . 

Figure [TD] plots structural parameters and [Fe/H] for our sample as a function of galac- 
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tocentric distance R gc . For the Milky Way clusters the true three-dimensional distance is 
used, while for external galaxies only projected distances are available. For the M31 clusters, 
no correlation of concentration c with R gc is seen. Central surface brightness fiyfi becomes 
fainter with increasing R gc but this is likely to be a selectio n effect: fainter clusters a re diffi- 
cult to see against the bright background of the M31 bulge ([McLaughlin et al.ll2007bl discuss 
the same effect in detail for NGC 5128). As expected, Rh increases with R gc , although the 
trend is largely driven by the innermost clusters in M31. In the Milky Way the Rh — R gc 
trend is driven both by small inner clusters and the large clusters at very large distances 
from the Galactic center. Such distant, low-surface brightness clusters will be difficult to 
detect in external galaxies, and their absence from the present sample may well be due to 
selection effects. For the c lusters in NGC 5128, which cover a more restricted range in R gc , 
McLaughlin et al.l (l2007bl ) conclude that the Rh — R gc trend is weak and only marginally 
significant. For the M31 and Milky Way clusters, a weak correlation of [Fe/H] with R gc is 
present, but again is driven mainly by clusters at the extremes of R gc . The Magellanic Cloud 
and Fornax dSph clusters are offset from the larger galaxies in all of these plots, as might 
be expected from th e parent galaxies' s maller sizes. Scaling cluster properties by galacto- 
centric distance as in iMcLaughlinl (120001 ) is therefore expected to be useful in computing the 
'fundamental plane' for these clusters. 

There are half a dozen very small (Rh < 1 pc) M31 clusters visible in the upper right 
panel of Figure [101 Such clusters are not found in the other galaxies studied, which makes 
their presence in M31 slightl y suspicious. Of th e six clusters, only the largest, B155, has a 
confirmatory radial velocity (IGalleti et al.ll2006l ) while the remainder are classified by those 
authors as "cluster candidates". The compactness of these objects makes assessing their 
nature problematic, even with the use of HST imaging. They are a reminder that the 
structural analysis of star clusters in distant galaxies will always be problematic: existing 
catalogs are neither fully complete nor entirely reliable. 

Figure [11] plots cluster structural parameters as a function of [Fe/H]; only clusters with 
measured values of [Fe/H] (not assumed average values) are shown. No correlation of concen- 
tration c with [Fe/H] is seen. A weak correlation of Rh with [Fe/H] is present in the expected 
sense: Rh decreases with increased metallicity. The exception is NGC 5128, for which the 
best-fit slope of log-R^ against [Fe/H] is zero; the present sample of clusters in this galaxy 
also has a lower metallicity gradient ([Fe/H] vs. R gc ) than in other galaxies. The Rh — [Fe/H] 
corrlation weakens if only clusters with 2 < R gc < 20 kpc are considered, suggesting that a 
different mix of bulge and halo clusters in M31 and the Milky Way compared to NGC 5128 
may be responsible. For central surface brightness /iy j0 there is a slight systematic increase 
with [Fe/H], likely because of the weak correlation of [Fe/H] and R gc . Finally, the lower-right 
panel shows that there is a weak mass-metallicity correlation: this is due entirely to our use 



of population sy nthesis model mass-to- ligh t ratios, which increase rapidly above [Fe/H] = — 1 
(see Figure 1 of McLaughlin et al.ll2007bl ). All of these trends, however, are mild compared 
with the correlations with cluster luminosity to be presented below. 



Figure [T2l shows some of the correlations between cluster parameters and luminosity, for 
GCs in M31, the Milky Way, NGC 5128, the Magellanic Clouds (MCs), and the Fornax dSph. 
The properties of clusters in all six galaxies fall in the same regions of parameter space. The 
selection effects discussed above are apparent in that there are fewer faint and low-surface- 
brightness clusters in the extragalactic samples. The new M31 ACS/STIS observations 
in particular do not contain any faint, low-concentration clusters, partly because the ACS 
observations targeted luminous M31 clusters and also because very low-concentration M31 
clusters are too well-resolved to be analyzed with the ellipse- fitting method used (see § 13. ip . 
There are few very large (Rh > 10 pc) extragalactic clusters, probably because these are 
difficult to recognize and more easily confused with background galaxies. The right-hand 
panels of the Figure show the lack of correlation between Rh and L and the strong correlation 
of Ef, with L which are both hallmarks of the globular cluster fundamental plane. The 
fundamental plane correlations evidently hold over nearly 4 decades in luminosity and in a 
variety of environments. 



6.3. The Fundamental Plane 



Correlations between globular cluster structural parameters are often summarized as 
describing a "fundamental plane" analogous to (but different from) that for elliptical galax- 
ies. There are several equivalent fo r mulations of the globular cluster fundamental plane in 



the lit erature (IBurstein et al.l 119971 ; iDjorgovskil 1 19951 ; McLaughlin! |2000| ; McLaughlin et al. 



2007bl ). Here we concentrate on the relationships between half-light radius, mass, and bind- 
ing energy to search for differences between the fundamental planes of clusters in different 
galaxies. An important caveat to this analysis is that the GC masses used are derived 
from population synthesis model-predicted mass-to-light ratios, not from directly observed 
velocity dispersions (see § 15. 3p . 

One of the many ways of looking at the globular cluster fundamental plane is a plot of 
projected cluster half-light radii as a function of mass. Modulo some stretching, this plot is 
equivalent to one of projected half-light surface density versus M, which in turn is the 
largest part of the K m) \ and parameter space. A different reformulation makes Rh vs. M 
equivalent to binding e nergy Eh vs. M: E h = f(c)M 2 /Rh where f(c) is a weak function of the 
cluster concentration (jMcLaughlinl |2000L and see below). Figure [TBI shows Rh versus M for 
GCs in all of our galaxies. As is well known from previous investigations of the Galactic GC 
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system (Ivan den Bergh et al.lll991t iDjorgovski fc Meylanlll994l ). there is no clear, monotonic 
mass-radius relation for globular clusters. However, at masses M > 1.5 x 10 6 t here is an 
apparent upturn in the lower envelope of Rh versus M. iMcLaughlin et al.l (l2007bl ) interpret 
this increasing lower bound on Rh for clusters as an extension of a very similar relation for 
early-type galaxies. The latter relati on follows from the ex istence of a "zone of exclusion" 
in the (luminosity-based) k space of iBurstein et al.l (119971 ) . which can be interpreted as a 
mass-dependent upper limit on the average surface density £/,. The bound i?™ m (M) = 
[M/27r£™ ax ] 1//2 for galaxies, simply extrapolated down into t he GC mass range, is drawn as 
the dash-dot line in Figure [T3| see (IMcLaughlin et al.ll2007bl ) for further details. 



The two dashed lines in the lower-left corner of the Figure represent a different bound 
on GC size: evaporation by two-bod y relaxation, which is exp ected to occur within 20-40 
initial half- mass relaxation times (e.g. lGnedin fc Ostrikerlll9971 ). The two lines correspond to 
t r h = (13/20) Gyr (upper) and t r h = (13/40) Gyr (lower); clusters which formed with masses 
and sizes to the lower-left of these lines are expected to have evaporated by the present day. 
Several of the M31 clusters are below the 'evaporation' lines, indicating that they must have 
evolved into this region of the diagram and may in fact be dissolving. The most extreme 
of these are the very compact M31 clusters discussed in § 16.21 The identfication of these 
objects as M31 GCs is somewhat suspect: high-resolution spectroscopic observations would 
be of interest to confirm their status and measure velocity dispersions. 

None of the galaxies' clusters show evidence for correlation of half-light radius with 
mass, but there are differences between galaxies in both position and scatter of their clusters 
in Figure [131 The larger scatter for Milky Way clusters compared to those in M31 and 
NGC 5128 can be attributed to the presence of more very diffuse c lusters in the MW sample 
(such clusters are difficult to detect in the extragalactic systems; IMcLaughlin et al.l l2007bl ) 
and greater distance uncertainties for individual Galactic GCs. The M31 sample has very 
few diffuse clusters due to the selection effects described above, and so is offset to smaller 
average Rh than the Milky Way or NGC 5128 samples The Magellanic Cloud and Fornax 
cluster samples are complete, however, and have larger average Rh than the large-galaxy 
GCs. 

Is the difference between cluster sizes in different galaxies due to local environment? The 
structure of tidally-limited clusters should depend on the local galaxy density, and cluster 
sizes are known to be correlated with galactocentric position (see the upper right panel of 
Figure [TO]) . A better indicator of galactic potential than galactocentric distance alone is: 



R 



Rgc/(8 k P c) 

V c /(220 km s" 1 ) 



(11) 
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where V c is th e galaxy circular velo city^] We 



use va 



220; M31, 230 dCarignan et al.l l2006h: NGC 5128. 190 dWqodlevl l20o3: IPeng et al. 



20041) . and Fornax dSph, 18 flWalker et al. 



ues V r . in km s 1 of: Milky Wa y. 



12004a 



Mathieu et al.l Il996l ): LMC. 65 (Ivan der Marel et all hood): SMC, 60 flStanimirovic et al. 



20061 ). For a galaxy whose total mass distribution 



is an isothermal sphere, mass density p oc (V c /R gc ) 2 , so R* c oc p~ l l 2 is a guide to the local 
galaxy density. Figure \\M shows Rh as a function of re-normalized galactocentric distance R* c , 
with the four galaxies plotted separately. The normalization has removed the Rh/R gc offset 
between the Magellanic Cloud/Fornax clusters and those in the large galaxies, indicating 
that the larger cluster sizes in the small galaxies can be understood as due to the lower 
galactic densities. The least-squares fit of log Rh = C + 7log(-R* c ) is shown for each galaxy, 
and given in Table [KB As might be expected, there is less dependence of Rh on i?* c for 
the extragalactic clusters compared to the Milky Way GCs, due to the use of projected 
galactocentric distances in computing R* c . 

The normalization by galactic potential should allow us to fa irly compare the fu nda- 
mental planes of c lusters in the different galaxies. In analogy with iMcLaughlinl (120001 ) and 
Harris et al.l (120021 ) we define the re-normalized quantities 



R*h = Rh x (R 



(12) 



and 



K 



E b x (R* gc y (13) 

where 7 is derived from the fits given above. Figure [15] shows that there is no correlation 
of R* h with [Fe/H], so any differences between galaxies are not due to residual effects of 
metallicity on the fundamental plane. 

The top panel of Figure [TBI shows El versus M, with the the corresponding least-squares 
fits for each galaxy given in Table [TO]. The basic proportionality Ef, oc M 2 is built into the 
definition of Eb, and clusters in all galaxies follow this relation very tightly. In particular, 
the fits for NGC 5128 and Milky Way clusters are identical within their errors, confirming 
the lHarris et al.l (120021 ) result with a sample three times larger and extending to much lower 
luminosities. At the high end of the mass range, the clusters in both M31 and NGC 5128 
fall below the best-fit lines, i ndicati ve of the 'break' in cluster properties at ~ 1.5 x 10 6 M Q 
noted by iMcLaughlin et al.l (j2007bl ). Below this break, the fits for all galaxies are similar, 
with the largest differences at the low-mass end where the number of clusters is small and 



The normalization of Rl r ensures tha t quantities derived from it are equivalent to those plotted for 



Milky Way clusters in IMcLaughlinl <|2000h . It should properly be computed with the three-dimensional 
galactocentric distance, but only the projected distance is available for the extragalactic clusters. 
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measurement errors are large. The bottom panel of Figure [16] shows the difference between 
El predicted from the Milky Way fit and the measured values. Table [10] gives the mean 
offsets for each galaxy. The NGC 5128 and Magellanic Cloud clusters (and the Milky Way 
clusters, by construction) have offsets consistent with zero, while the M31 clusters have a 
slight offset which can be attributed to sampling issues: the lack of diffuse clusters and the 
presence of the low-i?^ objects at M ~ 5 x 10 4 M Q 

Globular clusters have a strikingly similar nature in the widely different environments 
we have sampled, including a dwarf spheroidal, irregulars, large disk galaxies, and a giant 
elliptical. The homogeneous and internally precise data we can now work with, in addition 
to the much larger sample sizes of clusters, have established the trend of the fundamental 
plane securely, showing that globular clusters in all these environments follow the first-order 
trend Eb ~ M 2 quite accurately. The scatter around this basic relation is surprisingly small, 
showing that the structural properties of these star clusters are far simpler than we could 
have expected from theoretical arguments alone. Even the most massive GCs seem to show 
a fairly smooth transition between the properties of clusters and those of dwarf galaxies. 



7. Summary 



New Hubble Space Telescope observations are used to construct surface brightness pro- 
files for 34 M31 globular clusters, including the most massive clusters in the galaxy. Struc- 
tural models are fit to the surface brightness profiles: for reasons not yet understo od, M31 
clusters are unlike those in othe r near by galaxies in slightly preferring iKingi (119661 ) models 
over the more extended IWilsonl (119751 ) models. When combined with previously-published 
data, we now have a comprehensive and publicly available compilation of a complete suite 
of structural and dynamical parameters for as many GCs in M31 as in the Milky Way. The 
structural parameters of clusters in M31 and other local galaxies define a very tight fun- 
damental plane with similar slopes ~ M 2 , showing the essential similarity of clusters in 
different environments over a range of almost 4 decades in luminosity or mass. M31 clusters 
are slightly offset from the Eb — M relation defined by the Milky Way and NGC 5128, likely 
because of sampling effects. In both M31 and NGC 5128 clusters, the most massive GCs 
trend towards dwarf galaxies in the binding-energy/mass relation at about 10 6 M . The de- 
tailed modeling made possible by high-resolution imaging shows that the overall properties of 
globular clusters in different galaxies are remarkably similar; the subtle differences noted here 
may point to differences in the histories or environments of their parent galaxies. The present 
datas et pushes our assessment of the fundamental plane beyond the previous state of the 
field dMcLauehlinlboOol : Iflarris et al]l2002l : barmbv et aDbooi iMcLauehlin & van der Marel 



20051 ). making it possible to search for second-order trends in the FP line and probe the tran- 
sition between the structures of globular clusters and galaxies. 
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Fi g. 1— ACS/WFC F606W images of two M31 clusters. Left: B147, (classified as a star 
by lDubath fe Grillmairlll997l . based on a velocity dispersion measurement). Right: the very 
elliptical (e = 0.28) cluster B088. Both images have north up and east to the left and are 
square. 
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Fig. 2.— Histogram of population-synthesis ^-band mass-to- light ratios for M31 
GCs assuming an age of 13 G yr; and the same for Galactic globulars, from 



McLaughlin fc van der Marell (120051 ) . 
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Fig. 3. — Surface brightness profiles and model fits to four M31 GCs. Clusters are in order of 
increasing luminosity, from top to bot tom: B056D , B196, B023, a nd B 151. The thr ee panels 
for each cluster are, from left to right: iKingi (119661 ). IWilsonl (119751 ) and ISersid (119681 ) models. 
Dashed curves trace the PSF intensity profiles and solid (red) curves the PSF-convolved 
best-fit models with added background. Vertical dotted lines mark the radius where the 
best-fit cluster intensity is equal to the background. Open squares are ELLIPSE data points 
included in the least-squares model fitting and the asterisks are points (flags of BAD, SAT, 
or DEP in Table [2]) not used to constrain the fits. 
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Fig. 4. — Top: compa rison of structura l - mod el predictions for central velo city dispersion 
with o bservations from iDjorgovski et al.l (119971 ). iDubath fc Grillmairl (119971 ). and [Peterson 
(119891 ) . Open symbol is cluster B037, discussed in §5.51 Points connected by short-dashed 
lines are duplicate observations of the same cluster. Long-dashed line is line of equality; 
dot-dashed line is median ratio o" bs/o" P red = 0.85. Bottom: comparison of dynamical and 
population-synthesis derived mass-to-light ratios. Long-dashed line is line of equality; dot- 
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Fig. 5. — Comparison of parameters for model fits to 15 M31 clusters observed in both 
F606W and F814W filters. Top left: model total magnitude, top right: central surface 
brightness, bottom left: projected half-light radius, an d bottom rig ht: ratio of halfdi ght to 
core r adii (a measure of conc e ntrat ion) . Filled circles: iKingl (119661 ) models, squares: 
()1968l ) models, stars: IWilsonl (119751 ) models. 
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Fig. 7.— Comparison of IWilsonl (119751 ) and iKind (119661 ) model-fit parameters for M31 
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Errorbars show typical parameter uncertainties in a single model fit. 
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Fig. 8. — Top: distribution of ellipticity for globular clusters in the Milky Way (solid line), 
M31 (short-dashed lines), and NGC 5128 (long-dashed lines). Bottom: ellipticity versus 
projected galactocentric radius for globular clusters in M31 (filled squares), NGC 5128 (open 
circles), and the Milky Way (stars; three-dimensional R gc used). 
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Fig. 9. — Ellipticity versus luminosity and half-mass relaxation time for globular clusters in 
M31 (filled squares), NGC 5128 (open circles), and the Milky Way (stars). 
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Fig. 10. — Structural parameters as a function of galactocentric distance R gc for globular 
clusters in the Milky Way (stars), the Magellanic Clou ds and Fornax dwarf sphe roidal (open 
triangles), NG C 5128 [open circle s are clusters from McLaughlin et al.l l2007al ; small filled 
circles are from lHarris et al.l 120021 ] . and M31 (filled squares). The darker filled squares are 
the M31 clusters observed with ACS or STIS; lighter squares are WFPC2 observations. R gc 
is three-dimensional for Milky Way clusters but projected two-dimensional for other galaxies. 
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Fig. 11. — Structural parameters as a function of metallicity for globular clusters. Symbols 
as in Figure [TUl Only clusters with measured values of [Fe/H] (not assumed average values) 
are shown. 
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Fig. 12. — Structural parameters as a function of luminosity for globular clusters. Symbols 
as in Figure [TU1 
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Fig. 13. — Projected half-light radius as a function of cluster mass for GCs in the Milky 
Way, Magellanic Clouds and Fornax, M31, and NGC 5128. Symbols: Milky Way (stars), the 
Magellanic Cloud s and Fornax dwarf sphe roidal (open triangles), NG C 5128 [open circle s 
are clusters from McLaughlin et al.l l2007al : small filled circles are from lHarris et al.l 120021 ] , 
and M31 (filled squares). The darker filled squares are the M31 clusters observed with ACS 
or STIS; lighter squares are WFPC2 observ ations. Long-dashed lines are 'survival' lines; 
dot-dashed line is lower Rh bound defined by McLaughlin et al.l (j2007bl ). 
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Fig. 14. — Projected half-light radius as a function of re-normalized galactocentric distance 
R* c for globular clusters in the Milky Way (lower left), Magellanic Clouds and Fornax (lower 
right), M31 (upper right), and NGC 5128 (upper left). Solid lines are least-squares fits; 
dashed line in each panel is the Milky Way least-squares fit. 
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Fig. 15. — Re-normalized half-light radius R* h as a function of metallicity for globular clusters 
in the Milky Way, Magellanic Clouds and Fornax, M31, and NGC 5128. Symbols as in 
Figure [131 Only clusters with measured values of [Fe/H] (not assumed average values) are 
shown. 



-47- 




Fig. 16. — Re-normalized binding energy E£ (top) and difference from Milky Way least- 
squares fit (bottom), as a function of mass for globular clusters in the Milky Way, M31, 
and NGC 5128. Symbols as in Figure fT3| lines in top panel are least-squares fits for Milky 
Way (black solid line), all M31 clusters (grey solid line), M31 ACS/STIS clusters only (grey 
dashed line), NGC 5128 clusters (black dashed line), and Magellanic Cloud and Fornax 
clusters (dotted line). 
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Table 1. Integrated measurements for M31 clusters 



Cluster 


e 

[ 


6 

°E of N] 


V or alt. a F&lAW 
[VEGAMAG] [VEGAMAG] 


E(B-V) 


i?gc 
[kpc] 


[Fe/H] 


source b 


ACS/WFC targets from prof 


;ram 10260 












B037 


0.27 


51 


16.06 


14.07 


1.38 


2.97 


-1.07 ±0.2 


P02 


B042 


0.10 


134 


15.71 


14.29 


0.77 


3.31 


-0.78 ±0.2 


P02 


B063 


0.10 


28 


15.31 


14.13 


0.44 


3.47 


-0.87 ±0.3 


HBK 


B082 


0.05 


161 


15.02 


13.63 


0.72 


3.61 


-0.80 ± 0.2 


P02 


B088 


0.28 


178 


15.04 


13.90 


0.48 


3.79 


-1.81 ± 0.1 


P02 


B151 


0.07 


8 


13.26 


13.26 


0.38 


1.61 


-0.75 ±0.2 


HBK 


Serendipii 


ous clusters observed in program 10260 












B056D 


0.22 


117 


18.81 


17.72 


0.08 


4.20 


-1.20 ±0.6 


AVE 


B041 


0.05 


30 


18.95 


18.13 


0.13 


2.72 


-1.22 ± 0.2 


P02 


B061 


0.14 


2 


16.83 


15.62 


0.37 


3.58 


-0.73 ±0.3 


P02 


B090 


0.14 


124 


18.49 


17.52 


0.08 


3.15 


-1.20 ± 0.6 


AVE 


B147 


0.10 


172 


16.03 


15.04 


0.08 


1.29 


-0.24 ±0.4 


HBK 


B162 c 


0.30 


47 


17.62 


16.54 


0.31 


2.27 


-1.22 ± 0.2 


BOO 


ACS/HRC targets 


from programs 9719 and 9767 












B023 


0.13 


60 


13.88 


12.72 


0.36 


4.45 


-0.92 ± 0.1 


HBK 


B158 


0.05 


13 


14.36 


13.44 


0.12 


2.38 


-1.02 ±0.0 


P02 


B225 


0.05 


3 


13.87 


12.89 


0.15 


4.65 


-0.67 ±0.1 


P02 


G001 


0.20 


162 


13.79 




0.08 


34.55 


-0.95 ±0.1 


MEY 


STIS targets from 


program 8640 












B020 


0.05 


10 


14.97 




0.08 


7.37 


-1.07 ±0.1 


HBK 


B023 


0.12 


81 


14.12 




0.36 


4.45 


-0.92 ±0.1 


HBK 


B196 


0.05 


60 


17.36 




0.24 


8.10 


-2.46 ±0.6 


BOO 


B236 


0.16 


149 


17.44 




0.08 


8.49 


-1.20 ±0.6 


AVE 


B289 


0.13 


39 


16.27 




0.08 


16.88 


-1.20 ± 0.6 


AVE 


B336 


0.00 




17.91 




0.08 


12.91 


-1.20 ±0.6 


AVE 


B351 


0.05 


110 


17.61 




0.11 


12.57 


-1.43 ±0.2 


BOO 


B361 


0.16 





17.08 




0.08 


14.43 


-1.20 ±0.6 


AVE 


B363 


0.15 


173 


17.91 




0.08 


10.44 


-1.63 ± 0.4 


BOO 


B373 


0.09 


22 


15.63 




0.10 


10.07 


-0.50 ±0.2 


HBK 


B396 c 


0.18 


152 


17.45 




0.08 


17.27 


-1.20 ±0.6 


AVE 


B405 


0.09 


133 


15.25 




0.12 


18.20 


-1.80 ±0.3 


HBK 


B407 


0.06 


173 


16.11 




0.16 


19.76 


-0.85 ±0.3 


HBK 


B422 


0.41 


149 


17.98 




0.08 


16.34 


-1.20 ±0.6 


AVE 


B461 


0.16 


105 


17.44 




0.47 


10.87 


-1.36 ±0.5 


BOO 


B462 


0.09 


70 


18.03 




0.08 


3.74 


-1.20 ±0.6 


AVE 


G002 


0.09 


106 


15.93 




0.08 


33.62 


-1.70 ±0.4 


HBK 


G339 


0.10 


174 


17.35 




0.08 


28.68 


-1.20 ±0.6 


AVE 


G353 


0.16 


121 


17.31 




0.08 


26.32 


-1.20 ±0.6 


AVE 


WFPC2-observed clusters from Barmby et al. (2002) 












B006 


0.08 


74 


15.46 




0.13 


6.39 


-0.58 ±0.12 


HBK 


B011 


0.09 


76 


16.58 




0.08 


7.68 


-1.54 ±0.34 


HBK 


B012 


0.08 


16 


15.04 




0.11 


5.74 


-1.70 ±0.22 


P02 


B018 


0.15 


4 


17.53 




0.08 


9.27 


-1.63 ±0.77 


P02 


B027 


0.07 


97 


15.56 




0.16 


5.99 


-1.64 ±0.32 


HBK 
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Table 1 — Continued 



Cluster 


€ 


a 
u 

1 Hi Ol IN J 


v or dir. 
TVFP AM API 

1 V IZjV_xiT.lVliT.V_xl 


TVPP AM Afl 


njyiD — V j 


-Kgc 
[KpCJ 




source 


B030 


0.10 


118 


17.38 




0.66 


5.64 


-0.39 ± 0.36 


P02 


B045 


0.08 


40 


15.78 




0.18 


4.87 


-1.00 ± 0.26 


P02 


B058 


0.10 


138 


14.93 




0.12 


6.94 


-1.45 ± 0.24 


HBK 


B068 


0.22 


42 


16.41 




0.45 


4.31 


-0.29 ±0.59 


HBK 


B076 


0.09 


69 


16.93 




0.21 


2.85 


-0.73 ± 0.07 


P02 


B109 


0.10 


72 


16.20 




0.08 


1.39 


-0.17 ±0.49 


P02 


B110 


0.05 


49 


15.36 




0.16 


2.92 


-1.06 ± 0.14 


P02 


B114 


0.06 


132 


17.44 




0.16 


0.88 


-1.20 ± 0.60 


AVE 


B115 


0.08 


63 


16.00 




0.08 


0.64 


-0.15 ±0.38 


HBK 


B123 


0.14 


62 


17.42 




0.08 


1.28 


-1.20 ±0.60 


AVE 


B124 


0.07 


164 


14.78 




0.08 


0.21 


-1.20 ±0.60 


AVE 


B127 


0.08 


64 


14.47 




0.15 


0.32 


-0.87 ±0.17 


P02 


B143 


0.05 


158 


15.95 




0.13 


0.97 


0.09 ±0.42 


HBK 


B155 


0.12 


80 


18.01 




0.16 


3.14 


-1.20 ±0.60 


AVE 


B156 


0.05 


67 


16.97 




0.08 


3.60 


-1.51 ± 0.38 


P02 


B160 


0.18 


2 


18.08 




0.08 


3.56 


-1.17 ± 1.25 


P02 


B231 


0.17 


136 


17.25 




0.08 


5.54 


-1.49 ±0.41 


P02 


B232 


0.18 


42 


15.65 




0.07 


4.96 


-1.81 ±0.17 


P02 


B233 


0.11 


74 


15.72 




0.12 


8.05 


-1.59 ±0.32 


HBK 


B234 


0.07 


71 


16.78 




0.08 


6.00 


-0.95 ±0.13 


P02 


B240 


0.16 


98 


15.18 




0.08 


7.22 


-1.76 ±0.18 


HBK 


B264 


0.26 


142 


17.58 




0.08 


0.38 


-1.20 ±0.60 


AVE 


B268 


0.12 


103 


18.31 




0.08 


1.39 


-1.20 ± 0.60 


AVE 


B279 


0.20 


79 


18.55 




0.08 


7.74 


-1.20 ±0.60 


AVE 


B311 


0.09 


54 


15.45 




0.23 


13.06 


-1.93 ±0.08 


P02 


B315 


0.13 


159 


16.47 




0.08 


12.62 


-2.11 ± 0.53 


HBK 


B317 


0.11 


66 


16.57 




0.13 


10.14 


-2.12 ±0.36 


HBK 


B319 


0.0 





17.61 




0.08 


11.80 


-2.27 ±0.47 


BOO 


B324 


0.0 





18.45 




0.08 


8.29 


-1.20 ± 0.60 


AVE 


B328 


0.27 


159 


17.86 




0.08 


8.09 


-1.22 ±0.80 


HBK 


B330 


0.14 


102 


17.72 




0.08 


8.46 


-1.20 ±0.60 


AVE 


B331 


0.24 


70 


18.19 




0.08 


8.31 


-1.44 ±0.69 


BOO 


B333 


0.23 


26 


18.84 




0.08 


7.95 


-1.20 ±0.60 


AVE 


B338 


0.06 


102 


14.20 




0.12 


10.21 


-1.38 ±0.09 


HBK 


B343 


0.09 


70 


16.31 




0.07 


14.71 


-1.49 ±0.17 


HBK 


B358 


0.12 


63 


15.22 




0.08 


19.78 


-1.83 ±0.22 


HBK 


B368 


0.0 





17.92 




0.08 


9.52 


-1.20 ±0.60 


AVE 


B374 


0.21 


106 


18.31 




0.08 


9.64 


-1.90 ±0.67 


P02 


B379 


0.09 


55 


16.18 




0.10 


11.30 


-0.70 ±0.35 


HBK 


B384 


0.20 


121 


15.75 




0.12 


16.36 


-0.66 ±0.22 


HBK 


B386 


0.08 


140 


15.55 




0.08 


14.02 


-1.57 ±0.17 


P02 


B468 


0.0 





17.79 




0.08 


20.05 


-1.20 ±0.60 


AVE 


BH04 


0.32 


128 


19.69 




0.08 


8.33 


-1.20 ±0.60 


AVE 


BH05 


0.19 


123 


16.06 




0.08 


9.01 


-1.20 ±0.60 


AVE 


BH11 








19.86 




0.08 


9.41 


-1.20 ±0.60 


AVE 
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Table 1 — Continued 



Cluster 


€ 


e 


V or alt. a 


FS1AW 


E(B-V) 


i?gc 


[Fe/H] 


source* 3 






[°E of N] 


[VEGAMAG] 


[VEGAMAG] 




[kpc] 






BH18 


0.10 


116 


18.18 




0.08 


12.35 


-1.20 ±0.60 


AVE 


BH20 


0.18 


175 


18.07 




0.08 


4.44 


-1.20 ±0.60 


AVE 


BH21 


0.15 


30 


18.61 




0.08 


5.01 


-1.20 ±0.60 


AVE 


BH23 


0.09 


75 


18.83 




0.08 


1.28 


-1.20 ±0.60 


AVE 


BH24 


0.38 


14 


20.67 




0.08 


5.37 


-1.20 ±0.60 


AVE 


BH29 








19.65 




0.08 


8.71 


-1.20 ±0.60 


AVE 


DA038 


0.34 


108 


19.25 




0.08 


9.43 


-1.20 ±0.60 


AVE 


M091 








19.14 




0.08 


9.31 


-1.20 ±0.60 


AVE 


NB39 


0.13 


28 


17.94 




0.08 


0.19 


-1.20 ±0.60 


AVE 



a F606W for ACS-observed clusters, except G001 (F555W). 

b HBK= lHuchra et al ] jl99lh : P02= |Perrett et al ] (|2002h : B00= lBarmb7 et al AVE=no measurement; 

[Fe/H] = -1.2 ± 0.6 assigned. 

c These clusters appear close to the image edge; their measurements do not cover the full range of radial 
distance. 
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Table 2. 49 F606W, F814W, or V Intensity Profiles for 33 GCs in M31 



Name 


Detector 


Filter 


R 


/ 


uncertainty 


Flag 








r i 

[arcsec] 


r t — 9i 

[L Q pc 2 ] 


[L Q pc 2 ] 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


RD5RD 


WFC 

VV J. Vj 


FR1 4W 

J. OIt: VV 






1D1 ^84 

1U1 .OO'-r 


OK 


1 J\> U \) yJ 


WFC! 

VV J. 


F81 4W 

J. O 11 V V 




899 


1 07 fi99 


DFP 

yJ 1_jL 


B056D 


WFC 


F814W 


0.0315 


2233.302 


112.895 


DEP 


B056D 


WFC 


F814W 


0.0347 


2236.559 


117.452 


DEP 


B056D 


WFC 


F814W 


0.0381 


2241.757 


122.118 


DEP 


B056D 


WFC 


F814W 


0.0420 


2245.391 


127.393 


DEP 


B056D 


WFC 


F814W 


0.0461 


2245.646 


133.215 


DEP 


B056D 


WFC 


F814W 


0.0508 


2239.541 


138.504 


DEP 


B056D 


WFC 


F814W 


0.0558 


2234.688 


143.454 


OK 


B056D 


WFC 


F814W 


0.0614 


2229.836 


148.906 


DEP 


B056D 


WFC 


F814W 


0.0676 


2227.591 


154.365 


DEP 



Note. — Table [2] is published in its entirety in the electronic edition of 
the Journal. Only a small portion is shown here, for guidance regarding 
its form and content. See text for description of the FLAG column. 



Table 3. 


PSF Models: I PSF 




J [1 + (R/r ) a ] 


-0/a 


Detector 


Filter 


^0 


a 




FWHM 






[arcsec] 






[arcsec] 




WFC 


F606W 


0.0686 


3 


3.69 


0.125 


2.5 


WFC 


F814W 


0.0783 


3 


3.56 


0.145 


2.9 


HRC 


F555W 


0.0267 


3 


3.09 


0.0527 


2.1 


HRC 


F606W 


0.0294 


3 


2.95 


0.0593 


2.4 


HRC 


F814W 


0.0351 


3 


2.96 


0.0706 


2.8 


STIS 


CL 


0.074 


2 


4.2 


0.093 


1.8 



Table 4. Basic Parameters of Fits to 49 Profiles of 33 GCs in M31 



Name 


Detector 


Ext. 


V color 


Wpts 


Model 


\ 2 ■ 

' v min 


-fbkg 


Ho 


c/n 


Mo 


log r 


log r 






[mag] 


[mag] 








[L 0P c- 2 ] 






[mag arcsec - 2 ] 


[arcsec] 


[PC] 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 



B056D WFC/F606 0.224 0.274 ± 0.051 52 K66 

W 
S 

B056D WFC/F814 0.144 1.286 ± 0.064 52 K66 

W 
S 

B090 WFC/F814 0.144 1.114 ±0.064 50 K66 

W 
S 



Note. — See text for column descriptions. Table|4]is available in its entirety in the electronic edition of the Journal. A short extract from it is shown here, for guidance 
regarding its form and content. 



149.50 
131.14 
195.14 
175.21 
149.55 
224.33 
105.21 
107.84 
124.24 



20.38 
19.98 
20.92 
27.74 
27.25 
28.56 
76.95 
75.18 
77.42 



±0.76 
±0.97 
±0.60 
± 1.46 
± 1.55 
± 1.00 
±2.49 
±4.06 
±2.29 



6.20 
6.10 

5.80 
5.50 

7.50 
7.50 



+0.60 
0.60 
+0.70 
0.70 



+0.80 
0.90 
+ 0.80 



+0.60 
0.50 
+0.70 
0.70 



1.31 
1.84 
1.26 
1.21 
1.59 
1.14 
1 

2.79 
2.30 



+0.16 

0.15 
+0.40 
0.28 
+0.24 
0.20 
+0.21 
0.20 
+0.35 
0.33 
+0.28 
0.20 
+ 0.18 
-0.15 
+0.45 
0.55 
+ 0.60 
0.44 



17 nn+o.05 

17 05+0.05 

-+0.23 
-0.28 
3+0.08 
'-0.08 
,+0.07 
-0.06 
,+0.25 
-0.35 
,+0.09 
-0.11 
,+0.09 



17.45] 
16.1 
16.91] 
16.5lj 
15.71] 
15.76^ 
14.10 



0.08 
+ 0.60 



-0 562+ - 042 

U.ODZ_ Q Q42 

-0 526+ - 044 
U.OZD_ 043 

-0.896t°;^ 

-0 51 K+0.075 

u.oro_o 053 

-0.825±»:5S 
-0-898±g : g« 
-0.860t°;°g 

9 1QK+0-465 



016+ - 042 
U.U1D_ 042 

o n5i+ 044 

U.U01_ 043 
n o-|o+0.161 
u.aro_ 208 

016+ ' 068 

U.U1D_ 061 

063+ ' 075 
U.UDO_ 053 

-0.248t°; 2 6l 

-0 320+ 045 
u.ozu_ 05g 

_ n 909+0. 052 
u.zoz_ 053 

1 W +0.465 
1.00 I _ n f-Q A 



Table 5. Derived Structural and Photometric Parameters of GCs in M31 



Name Detector Model log r tid log R c log Rh \og(R h /R c ) log I log j log L v Vtot log//, 

[pc] [pc] [pc] [L 0i vpc~ 2 ] [L ,vpc -3 ] [Lq.v] [mag] [L Q ,vpc -2 ] 

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) 



B056D 


WFC/F606 


K66 


1 O 9 + 0.12 


-0. 


016 +0.032 
ulD -0.034 


343+ 057 





otro+0.092 
' o '- , °-0.070 


q 9 O+0.03 

°- zo -0.03 


2 99+ 06 
z,ya -0.06 


4 3 +0.04 
^• ou -0.03 


18 


r-,+0.08 
°^-0.10 


2 


00+O.O6 
° -0.08 






W 


1 89+ 036 
l.O3_ 24 





m9 +0.028 
uu -0.030 


n QCQ+o.102 

0.383_ 055 





ooi +0.132 
00 -0.083 


q 97 +0.03 
o - z '-0.03 


2 96+ 05 
z ' MD -0.05 


4 34+O.O6 
-0.04 


18 


45+"-io 
■^°-0.16 


2 


77 +0.08 
"-0.15 






S 


OO 


-0 


qio+0.161 
O1 °-0.208 


n q1 q+0.040 
U.O!O_ 031 





Kq1 +0.249 
00 -0.192 


q 47 + 0.11 

■^'-o.io 


o Qo+0.26 
^■ 8S -0.20 


4 26+ 03 
4 - ZD -0.03 


18 


oq+0.06 
"''-0.07 


2 


84+0-05 

°^-0.06 


B056D 


WFC/F814 


K66 


i 99 +0.15 
1 - zz -0.13 


-0 


n91 +0.049 
uz -0.049 


n 9 ao+0.058 
O.^!88_ 036 





o 10 +0.107 
Olu -0.085 


o nq+0.04 
°' za -0.04 


3 Q1+0.09 
°- u -0.09 


4 25+0' 04 
4 - zo -0.03 


18 


gg+0.08 
DD -0.10 


2 


oo+0. 06 
■88-0.09 






W 


1 66+ ' 29 

l.OD_ 25 





Q00 +0.042 
uuu -0.036 


u.ou»_ 045 





309+ ' 104 
oua -0.087 


q 9s +0.03 
■ 5 " ao -0.04 


9 07+0.07 
z - 3 ' -0.08 


4 28+ ' 05 
4 - z °-0.04 


18 


5O+0.10 

Jy -0.12 


2 


oc+0. 07 

■ 8b -o.io 






s 


OO 


-0 


948 +0.161 
z ^°-0.243 


266+"' 038 

U.ZDD_q 022 





dl -0.183 


3 44+0-14 


9 S9 +0.31 
z '° -0.20 


- 9 q+0.03 
^• zo -0.03 


18 


79 +0.07 
' ,z -0.07 


2 


g +0.04 
au -0.07 


B090 


WFC/F814 


K66 


1 36+° 13 


-0 


ooy+0.040 

•"'-o.oss 


259+ - 113 
v.z,oj_ u ()(iu 





cng+0.166 
Jyo -0.109 


q Oq+0.05 

■3-^-0.04 


q o fi +0.10 
3-»t>_ . 08 


4 40+0- 06 
4 -* u -0.04 


18 


■^-o.is 


3 


ns +0.10 






W 


9 c-i+0.40 


-0 


oqo+0.042 
oua -0.045 


445+0-329 
u -^°-0.198 





754+0-374 
■ ' J -0.240 


q oi +0.04 


q s9 +0.08 
o - oz -0.08 


4 51+0.15 
4 - 01 -o.io 


18 


n9 +0.26 
uz -0.36 


2 


o O +0.30 
• 8z -0.51 






s 


OO 


-1 


+0.465 
°°' -0.684 


259+ 0101 


1 


s1fi +0.785 
Sib -0.534 


4 4 7 +0.33 
-0.24 


4 oo+O- 92 
^ ,au -0.62 


4 38+0.06 


18 


qr+0.12 


3 


UD -0.15 



Note. — See text for column descriptions. Table [5] is available in its entirety in the electronic edition of the Journal. A short extract from it is shown here, for 
guidance regarding its form and content. 

co 



Table 6. Derived Dynamical Parameters of GCs in M31 



Name Detector T^° p Model log M to t log E b log S log p log S h log cr p .o log u csc ,o log t lh log f a 

[M Q L-\} [M q ] [erg] [M Q pc- 2 ] [M pc- 3 ] [M pc" 2 ] [km s" 1 ] [km s" 1 ] [yr] [M Q (pc km s" 1 )" 3 ] 

(1) (2) (3) ' (4) (5) (6) (7) (8) (9) (10) (11) (12) (13) 



B056D 


WFC/F606 


2 00 , 5 +0 - 617 
z.uuo_ 263 


K66 


4 fi0+ 012 


49.27t°;J 7 


, r; 8 +0.12 


3. 


9O+0.13 
ZM -0.08 


„ -,9+0.13 
c> - lz -0.10 


52fi +0 ° 59 

U.0ZD_q 033 


1 


1 01 +0-059 
lul -0.033 


,,,+0.12 


0.484l°;J° 3 








W 


4 64+ 13 


49.29l°;J 7 


o c-7+0.12 
°- o '-0.07 


3 


9K+0.13 
ZD -0.08 


3 07+° 14 


n ™+0.059 
U.OOU_ 033 


1 


llfi +0.059 
11D -0.033 


+0.19 
8-Ot-O.lO 


0.434l°;° 7 7 








s 


4 5fi+° 12 


49.24l°;J 7 


o 77 +0.16 
ci -' '-0.11 


3. 


ic+0.27 
lo -0.20 


o U +0.13 


n cin+o.oeo 


1 


-,90+0.059 
• lzJ -0 .032 


s , 7 +0.09 
o- ' -0.06 


0.433l°;^ 77 


B056D 


WFC/F814 


Z.UUO_ 263 


K66 


4 - OD -0.07 


49.24l°;£ 7 


3 R0+ ' 12 


3 


O-.+0.15 

■•"-en 


o -,0+0.13 


n t-90+0.060 
u.dza_ 034 


1 


nqc; +0.060 
u3,J -0.035 


01 +0.12 
°-'-' 1 -0.07 


0.485l°;« 7 








W 


4 ^S +0 ' 13 


49 25+ - 17 


3 58+ ' 12 


3. 


o7 +0.13 
z ' -0.10 


o i 7 +0.13 
-0.12 


n co 4 +0.060 

U.0J4_ Q 034 


1 


11n +0.060 
llu -0.034 


s oc+0.13 
< 5 " :>o -0.09 


0.423l°;i°! 








s 


^' OJ -0.07 


49.22t°;J 7 


o 74 +0.18 
J ' ' 4 -0.12 


3 


, 9 +0.31 
±z -0.20 


3 2n+ ' 12 

c> - zu -0.09 


0.5241°;°!? 


1 


,91+0.060 
lzl -0.034 


o7 +0.09 
8 - z '-0.05 


0.3241°;^ 


B090 


WFC/F814 


2 nfK +0 - 617 

z.uuo_q 263 


K66 


4 7n+° 13 

4 - ' u -0.08 


49.56±°; 7 


4 1Q+0.13 
^■ iJ -0.07 


1. 


lfi +0.15 
iD -0.10 


o OO+0.16 
o.oo_ lg 


0.63818;°!° 


1 


24f) +0.060 
z4u -0.034 


s Qo+0.20 
°- oo -0.12 


1 (,01 +0.132 
i.uor_ og6 








W 


4 ai +0-19 
4.81_ 12 


49.59t°; 1 1 7 


4 12+ 012 
4 - iz -0.07 


4. 


in+0.14 
1 -0.10 


o 19 + 0.32 

°- lz -0.52 


n fi 4Q+o.060 
U.O4J_ 034 


1 


9c;4 +0. 060 
• zo ^-0.034 


8.65t°;«S 


0.962±°;i 2 4 








S 




49.52±°;J 7 


A 7 o+0.35 
4 -< 8 -0.25 


5. 


9O+0.92 
zu -0.63 


3 36+° ' 15 
°- oo _0 .16 


n TO +0.074 


1 


974 +0.062 
• z '*-0.036 


09+0. 18 
"■•"-0.12 


^- ' u ' -0.929 



Note. — See text for column descriptions. Table [6] is available in its entirety in the electronic edition of the Journal. A short extract from it is shown here, for guidance Iregarding its 
form and content. ^ 



Table 7. Galactocentric Radii and K-Space Parameters for GCs in M31 



Name Detector R gc Model Km 1 Km 2 Km 3 

[kpc] 

(1) (2) (3) (4) (5) (6) (7) 



B056D 


WFC/F606 


4.20 


K66 


-1. 


1Qfi +0.089 
lou -0.049 


4 059+° 158 
t.uoa_g 124 


341+ 007 
"■^ -0.002 








W 


-1. 


1Q1 +0.105 
iUi -0.055 


4011 +0.169 
*' ull -0.183 


347+0 017 
u -°^' -0.005 








S 


-1. 


1R o+0.084 
J-08-0.047 


4.085±°;^ 


33 6 +0.00i 


B056D 


WFC/F814 


4.20 


K66 


-1. 


16g +0.089 


4 136+ ' 160 


340+0.006 
u.oiu_ 0()0 








W 


-1. 


147 +0.093 
±fl ' -0.051 


4 120+ ' 164 
<i.±zu_ 13g 


n 043+0.008 

U.d4d_ 0()1 








s 


-1. 


1Q9 +0.084 
13 -0.048 


4 l% +al52 

1.J.OD_q 113 


n oo 7 +0.000 

u.oo/ _ 005 


B090 


WFC/F814 


3.15 


K66 


-1. 


UJD -0.064 


4 401+ ' 185 
, '* u -0 .207 


366+ ' 027 

U.JDD_ q 14 








W 


-0. 


oQo+0.239 

»y»_ .i4i 


— 0. ■>(>■> 


415+0-100 








s 


-1. 


icq+0.085 
1Ja -0.055 


4 315+0.205 


n 977 +0.033 
u - z ' ' -0.054 



Note. — See text for column descri p tions, particularly difference between n definitions 
used here and those of iBender et al. I l ll992l) . Table [7] is available in its entirety in the 
electronic edition of the Journal. A short extract from it is shown here, for guidance 
regarding its form and content. 



Table 8. Predicted Aperture Velocity Dispersions for GCs in M31 



Name Detector TP op Model log R h log R h log a ap (R h /8) log <r ap (K h /4) log cr* p (R h ) log <r ap (4il fe ) log <7 ap (8 iJ h ) 

[M@ Lq v ] [pc] [arcsec] [km s _1 ] [km s~ x ] [kms^ 1 ] [km s" 1 ] [km s _1 ] 

(1) (2) (3)'' (4) (5) (6) (7) (8) (9) " (10) (11) 



B056D 


WFC/F606 


2 no5+°- 617 

z.uuo_ 263 


K66 


0.343 


n 9 nc-+0.057 
u.zoo_ 038 





c 9 k+0.059 
OZO -0.033 





c 9 o+0.059 
oz,;i -0.033 


5m +0-059 
U.OU1_ 033 





4 ,- 4 +0.060 
■ 4O4 -0.033 


n 44R+0-0 60 
U.44D_ 34 








w 


0.383 


-0 1Q4+ 102 
u - iy4 -0.055 





c-nQ+0.059 
OZM -0.033 





cofi+0.059 
OZD -0.033 


n 4OQ+0 060 

u.*aa_ o 3 4 





4 „+0.060 
4oc> -0.034 


n 441+0-0 60 
u - 441 -0.036 








s 


0.313 


-0 265+ ' 040 

U.ZDO_q q 31 





4q4 +0.060 
4y4 -0.033 





p-r.o + 0.060 

OU °-0.033 


n 5n o+0.060 
u.ouo_ 033 





4 ,- 4 +0.060 
' 4o4 -0.034 


n 45n+ -0 60 
U.40U_ 34 


B056D 


WFC/F814 


2 oo5+°- 617 

z.uuo_ 263 


K66 


0.288 


n 9SQ +0.058 

-u.2»y_ 036 





cno+0.060 
Oz °-0.034 





r97 +0.0(i() 
oz ' -0.035 


n 5f)5 +0.060 
u.ouo_ 035 





4C ro+0-060 
■ 4t>8 -0.035 


n 459 +0.061 
U,4JZ -0.036 








W 


0.309 


-0 269+ 068 
u.zoy_ 045 





coo+0.060 
o,3o -0.034, 





™ +0.060 
°'-' 1 -0.034 


0.506±°;g|g 





4fif) +0.060 
4DU -0.035 


n 450+ lu,fil 
U.4OU_ 036 








s 


0.266 


n oi, +0.038 
U.J11_ 22 





r m +0.061 
OU1 -0.036 





^+0.061 
ol4 -0.035 


0.5081°;°°° 





45 q+0-060 
409 -0.035 


n 45fi+o- 06 o 

U.4OD_ 3 6 


B090 


WFC/F814 


2 no5+°- 617 

z.uuo_ 26 3 


K66 


0.259 


n qi o+0. 113 
U.OlO_ 69 





fi o R +0.060 
D,; ' D -0.034 





fi on+0.060 
D '-'' ;i -0.034 


0.6081°°°° 





ct; fi +0.060 
OOD -0.037 


n E 42 +0.061 
u.otz_ 038 








W 


0.445 


n 1,9+0.329 
u.ijz_ lgg 





K o Q +0.060 
D '-' M -0.036 





fi o 9 +0.060 
O'-' -0.039 


n 505+O.O62 
u.ovo_ 047 





.,,0+0.063 
o4,;> -0.053 


n r 24 +0.06S 
U.OZ4_ 055 








s 


0.259 


n Ql R+0-101 
U.OlO_ 068 





fina+0-060 
bUS -0.033 





fi91 +0.060 
' DZ -0.033 


601+°' ()(i0 
U.DU1_ 035 





548+0-060 
o4 °-0.036 


n cq 7 +0.061 
U.JOI_o 037 



Note. — See text for column descriptions. Table [8] is available in its entirety in the electronic edition of the Journal. A short extract from it is shown here, 
for guidance regarding its form and content. 
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Table 9. Measured Aperture Velocity Dispersions for GCs in M31 



Name 


Tap 


Aperture 


Source 




[km s _1 ] 


[arcsec] 





Clusters in the present sample 



B006 


11.56 ±0.11 


1.60 


B006 


10.6 ±0.4 


1.80 


B012 


16.15 ±0.45 


1.60 


B020 


14.27 ±0.24 


1.60 


B020 


15.3 ±0.5 


1.55 


B020 


18 ±3 


1.2 


B023 


25.46 ±0.90 


1.60 


B023 


24 ±6 


1.2 


B037 


19.55 ±3.50 


2.77 


B045 


9.82 ±0.18 


1.60 


B045 


8.7 ±0.5 


1.80 


B082 


25 ±5 


1.2 


B147 


11 ±2 


1.2 


B158 


20.50 ±0.26 


1.60 


B158 


21.9 ± 1.3 


1.55 


B158 


40 ±6 


1.2 


B225 


25.94 ±0.38 


1.60 


B225 


26.9 ±0.5 


1.80 


B225 


< 40 


1.2 


B240 


11.92 ±0.22 


1.60 


B289 


8.43 ±0.44 


1.60 


B343 


9.08 ±0.39 


1.60 


B343 


10.2 ± 1.7 


1.55 


B358 


8.11 ±0.36 


1.60 


B358 


7.1 ± 1.8 


1.55 


B373 


12.58 ±0.15 


1.60 


B379 


8.15 ±0.20 


1.60 


B384 


10.10 ±0.16 


1.60 


B384 


9.1 ± 0.5 


1.55 


B386 


11.49 ±0.24 


1.60 


B405 


8.57 ±0.45 


1.60 


B407 


9.52 ±0.11 


1.60 


G001 


25.06 ±0.32 


1.60 


G001 


21.47 ± 1.60 


2.77 


G002 


9.70 ±0.29 


1.60 


Clusters not in the present sample 


B019 


19 ±3 


1.2 


B163 


21 ±4 


1.2 


B171 


18 ±3 


1.2 


B193 


13.20 ±0.24 


1.60 


B193 


12 ±3 


1.2 


B218 


17.62 ±0.24 


1.60 


B218 


16.3 ±0.8 


1.55 



Diorgovski et al. (1997) 
Dubath fe Grillmair (19971 
Diorgovski et al. (1997) 
Diorgovski et al. (1997) 
Dubath fe Grillmair (19971 
Peterson (1989) 
Diorgovski et al. (1997 1 
Peterson (1989) 



Cohen (20061 
Diorgovski et al. (1997 1 
Dubath fe Grillmair (19971 
Peterson (19891 
Peterson (19891 
Djorgovsk^et^aL (1997) 
Dubath fe Grillmair (19971 
Peterson (19891 
Diorgovski et al. (1997 1 
Dubath fc Grillmair (19971 



Peterson (19891 
Diorgovski et al. (1997 1 
Diorgovski et al. (1997 1 
Diorgovski et al. (1997 1 
Dubath fc Grillmair (19971 
Diorgovski et al. (1997 1 
Dubath fc Grillmair (1997) 
Diorgovski et al. (1997 1 
Diorgovski et al. (1997 1 
Diorgovski et al. (19971 
Dubath fc Grillmair (19971 
Djorgovsk^et^L 
Diorgovski et al. 
Diorgovski et al. 
Diorgovski et al. 
Cohen (20061 
Diorgovski et al. 



Peterson (1989) 
Peterson (19891 
Peterson (19891 
Diorgovski et al 
Peterson (1989 
Diorgovski et al 
Dubath fc Grillmair (19971 
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Table 9 — Continued 



Name 


Tap 


Aperture 


Source 




[km s- 1 ] 


[arcsec] 




B224 


< 10 


1.2 


Peterson f 19891 


B381 


< 10 


1.2 


Peterson (1989) 



Table 10. Fundamental plane fits 



Galaxy 


log R h = 


C + 7 log(i?* c ) 




= a(logM) + b 










C 


7 


RMS 


a 


b 


RMS 


5\ogE* 


N 


Milky Way 


0.57 ±0.02 


0.47 ±0.05 


0.21 


2.04 ±0.04 


39.63 ±0.21 


0.21 


0.02 ±0.02 


85 


NGC 5128 


0.58 ±0.02 


0.19 ±0.07 


0.19 


2.09 ±0.04 


39.36 ±0.20 


0.19 


-0.01 ±0.02 


104 


M31 


0.43 ±0.02 


0.20 ±0.04 


0.17 


1.92 ±0.03 


40.43 ±0.15 


0.17 


0.15 ±0.02 


84 


M31(ACS/STIS) 


0.47 ±0.02 


0.11 ±0.04 


0.09 


1.96 ±0.03 


40.19 ±0.15 


0.10 


0.09 ±0.02 


33 


MCs+Fornax 


0.65 ±0.04 


0.41 ±0.10 


0.17 


2.25 ±0.12 


38.48 ±0.65 


0.15 


-0.07 ±0.04 


18 



